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I. INTRODUCTION 

In the Standard Model (SM), semileptonic decays of 
B mesons proceed via first-order electroweak interactions 
and are mediated by the W boson [3-01. Decays involv- 
ing electrons and muons are expected to be insensitive to 
non-SM contributions and therefore have been the bases 
of the determination of the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix elements |V c f,| and \V u b\ [J]- Decays in- 
volving the higher-mass r lepton provide additional in- 
formation on SM processes and are sensitive to addi- 
tional amplitudes, such as those involving an interme- 
diate charged Higgs boson 0-0]. Thus, they offer an ex- 
cellent opportunity to search for this and other non-SM 
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contributions. 

Over the past two decades, the development of heavy- 
quark effective theory (HQET) and precise measurements 
of B -> D^tTVi decays [l(3| at the B factories [TJ HI 
have greatly improved our understanding of exclusive 
semileptonic decays. The relative rates 



K(D) 



B(B -» Dt-V t ) 
B(B ->• Dl-V t ) 



K(D* 



B(B -» D*t-V t ) 
B(B — s- D* 



(1) 

are independent of the CKM element \V c b\ and also, to 
a large extent, of the parameterization of the hadronic 
matrix elements. SM expectations @ for the ratios 
TZ(D) and 1Z(D*) have uncertainties of less than 6% and 
2%, respectively. Calculations based on two-Higgs- 
doublet models predict a substantial impact on the ratio 
1Z(D), and a smaller effect on 1Z(D*) due to the spin of 
the D* meson. 

The decay B D*t~V t was hrst observed in 2007 
by the Belle Collaboration [l3|. Since then, both BABAR 
and Belle have published improved measurements, and 
have found evidence for B Dt~~D t decays [143 - fl6j . Up 
to now, the measured values for TZ(D) and 1Z(D*) have 
consistently exceeded the SM expectations, though the 
significance of the excess is low due to the large statistical 
uncertainties. 
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We recently presented an update of the earlier mea- 
surement Q based on the full BABAR data sample fl7| . 
This update included improvements to the event recon- 
struction that increased the signal efficiency by more 
than a factor of 3. In the following, we describe the anal- 
ysis in greater detail, present the distributions of some 
important kinematic variables, and expand the interpre- 
tation of the results. 



W~/H , 




FIG. 1. Parton level diagram for B — > D^*'t~V t decays. 
The gluon lines illustrate the QCD interactions that affect 
the hadronic part of the amplitude. 



We choose to reconstruct only the purely leptonic de- 
cays of the t lepton, t~ — > e~z7 e j/ T and r~ — > pTV il y T , 
so that B -> D^t~V t and B -> D^£~Vi decays are 
identified by the same particles in the final state. This 
leads to the cancellation of various detection efficiencies 
and the reduction of related uncertainties on the ratios 
K(HW). 

Candidate events originating from Y(4S) —> BB de- 
cays are selected by fully reconstructing the hadronic de- 
cay of one of the B mesons (Stag), and identifying the 
semileptonic decay of the other B by a charm meson 
(charged or neutral D or D* meson), a charged lepton 
(either e or /i) and the missing momentum and energy in 
the whole event. 

Yields for the signal decays B — > D^t~V t and the 
normalization decays B — > Z)M £~Ve are extracted by an 
unbinncd maximum-likelihood fit to the two-dimensional 
distributions of the invariant mass of the undetected par- 
ticles m^ liss = p 2 niss = (p e +e-~PB tas ~PD(*)-Pe) 2 (where 
Pe+e-j P-Btagj and pi refer to the four-momenta of 

the colliding beams, the Stag, the D^*\ and the charged 
lepton, respectively) versus the lepton three-momentum 
in the B rest frame, \p* g \. The m^ iss distribution for de- 
cays with a single missing neutrino peaks at zero, whereas 
signal events, which have three missing neutrinos, have a 
broad w^iss distribution that extends to about 9GcV 2 . 
The observed lepton in signal events is a secondary par- 
ticle from the r decay, so its |pt | spectrum is softer than 
for primary leptons in normalization decays. 

The principal sources of background originate from BB 



decays and from continuum events, i.e. 



-> mi) 



pair production, where / = u,d,s,c,T. The yields and 
distributions of these two background sources are derived 
from selected data control samples. The background de- 
cays that are most difficult to separate from signal decays 
come from semileptonic decays to higher-mass, excited 
charm mesons, since they can produce similar and 
\Pi | values to signal decays and their branching fractions 
and decay properties are not well known. Thus, their 
impact on the signal yield is examined in detail. 



II. THEORY OF B 



DECAYS 



Standard Model 



Given that leptons are not affected by quantum chro- 
modynamic (QCD) interactions (see Fig. [T]), the matrix 
element of B — > D^t~V t decays can be factorized in 
the form [j| 

GfVci, 



(2) 

where L\ T and H, dM are the leptonic and hadronic 

Aw Aw * 

currents defined as 
L X ^ w (q 2 ,e T ) ee e M (A w )(rF T |r 7 /1 (l-75)^|0), (3) 
ff^<"V) = e;(X w )(D^\cY(l-7B)b\B). (4) 

Here, the indices A refer to the helicities of the W, D^*\ 
and r, q = pb — Pd(*"> i s the four-momentum of the virtual 
W, and T is the angle between the r and the three- 
momenta measured in the rest frame of the virtual W . 
The metric factor r\ in Eq. [5] is J?{±,o,s} = {1, 1, — 1}, 
where Xw = ±, 0, and s refer to the four helicity states 
of the virtual W boson (s is the scalar state which, of 
course, has helicity 0). 

The leptonic currents can be calculated analytically 
with the standard framework of clcctroweak interactions. 
In the rest frame of the virtual W (W*), they take the 
form [Hi: 



(5) 
(6) 
(7) 



L± = —2y // ~q 2 ~vd±, L± = ^fV2m T vdo, 



L = —2\/q*vdo 1 Lq = \/2m T v(d + — dJ), 
L~ = 0, = —2m T v, 



with 



/ nil I ± cos 9 T 

I - — , d± = — -j= — , do = sin0 r . (8) 



The choice of the selection criteria and fit configura- 
tion are based on samples of simulated and data events. 
To avoid bias in the determination of the signal yield, 
the signal region was blinded for data until the analysis 
procedure was settled. 



Given that the average q 2 in B — > D^t v t decays is 
about 8 GeV 2 , the fraction of t~ leptons with positive 
helicity is about 30% in the SM. 

Due to the nonperturbative nature of the QCD inter- 
action at this energy scale, the hadronic currents cannot 



7 



be calculated analytically. They arc expressed in terms 
of form factors (FF) as functions of q 2 (see Sees. Ill A II 
and|HA2J). 

The differential decay rate, integrated over angles, is 
derived from Eq. [2] and Eqs.EHZ] : 



dr T 



g% \v cb \ 



\p 



967T 3 TO^ 

l«oR(i + g 



3m? 



(I#4 



l-ff-l 



2q 2 



(9) 



where | is the three-momentum of the D^*> meson 

in the B rest frame. For simplicity, the helicities of the 
Z)W meson and the q 2 dependence of the hadron helicity 
amplitudes -ff±,o,s have been omitted. The assignment is 
unambiguous because in B — > D*t~V t decays, H± only 
receive contributions from Ad» = ±, while Hq s require 
\d* =0. In B — > Dt~V t decays, only = s is possible, 
which implies H± = 0. 



where ro» = /(m B + to_d*)- Using dispersion 

relations and analyticity constraints [3, fl9j , the universal 
FF and the ratios can be expressed in terms of just five 
parameters: 

h Al (w) = h Al (l) [1 - 8p 2 D ,z(w) + (53 P 2 D , - 15)z(io) 2 

-(231 P 2 D , -91)z(«;) 3 ], 
fli(io) = i?i(l) - 0.12(io - 1) + 0.05(io - l) 2 , 
R 2 (w) = R 2 (l) + 0.11(u; - 1) - 0.06(w; - l) 2 , 
R (w) = R (l) - 0.11(w; - 1) + 0.01(io - l) 2 . 

Here, z{w) = ( y/wTl- v / 2)/( \JwTT+ \/2). The factor 
h Al (1) only affects the overall normalization, so it cancels 
in the ratio K(D*). 

Three of the remaining four FF parameters, i?i(l), 
i?2(l), and p 2 D * , have been measured in analyses of 
B —¥ D*£~Vi decays. The most recent averages by the 
Heavy Quark Averaging Group (HFAG) Q and their cor- 
relations C are: 



1. Form factor parameterization of B —¥ D*t v t decays 

Four independent FFs, V, A , Ai, and A 2 , describe 
the non-perturbative QCD interactions in B — > D*t~v t 
decays. Based on the FF convention of Ref. @, the 
hadronic currents take the following form: 



H±(q 2 ) 

H (q 2 ) 



(m B +m D ,)A 1 {q 2 ) =F 



2m B 



1 



2m d* v q 



im%\p D ,\ 



m B + ra D * 



m B + mo 



■A 2 (q 2 ) 



■\P*D*\V(q 2 ), 



{m 2 B -m 2 D , - q 2 )(m B + m D ,)A 1 (q 2 ) 



H s (q 2 



2m B \p* D , 



■A (q 2 



(10) 



In this analysis, we use an HQET-based parameteriza- 
tions for the FFs that is expressed in terms of the scalar 
product of the B and D* four- velocities 



10 = v B ■ v D * 



2 2 

m D , — q" 



2mrj*m B 



(11) 



Its minimum value io m ; n = 1 corresponds to q^&x — 
(m B — mn*) ■ The maximum value is obtained for the 
lowest possible value of q 2 , which is the square of the 
mass of the lepton. Thus, io max = 1.35 for B -4 D*t~T> t 
decays and io max = 1.51 for B —> D*l~vi decays. 

In this framework, the FFs are usually expressed in 
terms of a universal form factor h Al (w) and ratios Ri (10): 

A\[w) = — - — r D -h Al (w), A {w) = h Al (w), 



r D * 



R 2 (w) 
A 2 (w) = n, Al (w), 



Ri(w) 
V(w) = h Al {w), 



p 2 D , =1.207 ±0.028, C{p 2 D ,,Ri{\)) =0.566, 
Rx{l) =1.401 ± 0.033, C(p 2 D ,,R 2 (l)) =- 0.807, 
R 2 (l) =0.854 ± 0.020, C(R X (1), R 2 (l)) =- 0.758. 

Ro(w) affects the decay rate only via the scalar 
hadronic amplitude H s (q 2 ). The corresponding leptonic 
amplitude L s (q 2 ,9 T ) is helicity suppressed, i.e., its rate 
is proportional to the mass of the lepton (Eq. [6]) . As a 
result, B —> D*£~vi decays are not sensitive to this FF, 
and Rq(w) has not been measured. We therefore rely 
on a theoretical estimate, Rq(1) = 1.14 ± 0.07, based on 
HQET @. 



2. Form factor parameterization of B — ^ Dr v T decays 

The non-perturbative QCD interactions in B — > 
Dt~v t decays are described by two independent FFs, 
referred to as V\ and S\ |||. The helicity amplitudes 
take the form: 



H (w) —^m B m D 



m B + m D 



y/qHw) 



a/io 2 - lVi(io), (12) 



H s (w) =yJm B m D mB mD (10 + 1) ffi(io). (13) 

The amplitudes corresponding to the helicities Xw = ± 
vanish because the D meson has spin 0. For this decay 
mode, the variable 10 is defined as in Eq. 111! except that 
the D* meson mass is replaced by the D meson mass mr>. 

Taking into account dispersion relations [l{|, V\ can 
be expressed as 

Vi(w) = Vi(l) x [1 - 8p 2 D z(w) + (51p 2 D - lO)z(io) 2 

-(252p 2 D -84)z(w)% (14) 
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where Vi(l) and p 2 D are FF parameters. The normal- 
ization Vi (1) cancels in the ratio TZ(D). Based on B — > 
Dl~Vn decays, the average value of the shape parameter 
is p 2 D = 1.186 ± 0.055 Hf. As for B -> D*t~V t decays, 
the scalar hadronic amplitude is helicity suppressed and 
as a result, S\{w) cannot be measured with B — > Di~Vi 
decays. We use instead the following estimate based on 
HQET [|: 



Sx(w) =V 1 (w){l + A[- 0.019 + 0.041 (w - 1) 
- 0.015(i<;- l) 2 ]}, 



(15) 



with A = 1 ± 1. 

We have employed this FF parameterization to gener- 
ate B — > Dt~V t and B — > Dt~Vi decays, as described 
in Sec. IIII C 21 Though we used the same FF definitions 
and parameters, we found a difference of 1% between the 
value of 1Z(D) that we obtained by integrating Eq. [S]and 
the value quoted in Ref. @. 

On the other hand, if we adopt the FF parameters of 
Ref. [2(|, we perfectly reproduce the TZ(D) predictions 
presented there. The translation of the FF parameteri- 
zation of Ref. (20| into standard hadronic amplitudes is 
not straightforward, so we do not use these FFs in the 
Monte Carlo simulation. Since both parameterizations 
yield essentially identical q 2 spectra, they are equivalent 
with respect to Monte Carlo generation, which is not sen- 
sitive to differences in normalization. 



3. SM calculation of IZ(D^) and q 2 spectrum 

We determine the SM predictions for the ratios 
TZ(D^*') integrating the expression for the differential de- 
cay rate (Eq. ^ as follows: 



B(B -> D^tv) _ S%2 
B(B -> £)(*)&/) ~ 7^1 

J mi 



dr T 



dq 2 



dr. 
dp 



dq 2 



(16) 



with ?max = ( m B - m D{ ,)) 2 . 

The uncertainty of this calculation is determined by 
generating one million random sets of values for all the 
FF parameters assuming Gaussian distributions for the 
uncertainties and including their correlations. We calcu- 
late 1Z(D(*') with each set of values, and assign the root 
mean square (RMS) of its distribution as the uncertainty. 
We apply this procedure for B° and B~ decays, and for 
£ = e and p, and average the four results to arrive at the 
following predictions, 



K(D) SM = 0.297 ±0.017, 
ll(D*) SM = 0.252 ±0.003. 



(17) 
(18) 



Additional uncertainties that have not been taken into 
account could contribute at the percent level. For in- 
stance, some electromagnetic corrections could affect 
B -> D^fvi and B -> D^t~V t decays differently 
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FIG. 2. (Color online). Predicted q 2 spectra for (a) B — > 
Dt~V t and B -¥ Dl'Vj decays for Vi(l)V cb = 0.0427 and 
(b) B ->• D*t~v t and B -¥ D*£~Vt decays for h Al {l)V cb = 
0.0359 il. 



The experimental uncertainty on TZ(D^) is expected to 
be considerably larger. 

The q 2 spectra for B — > LK*V _ F T decays in Fig. [2] 
clearly show the threshold at g^- = m 2 , while for 
B -> D(*H~Vt decays ~ 0. We take advantage 
of this difference in the signal selection by imposing 
q 2 > 4 GeV 2 . The spectra for I = e and p arc almost 

- 0.011 GeV 2 . 



identical, except for q 2 < m 2 



B. Two-Higgs-Doublet Model Type II 

As we noted in the introduction, B — > D(*'t~V t de- 
cays are potentially sensitive to new physics (NP) pro- 
cesses. Of particular interest is the two-Higgs-doublct 
model (2HDM) of type II, which describes the Higgs sec- 
tor of the Minimal Supersymmetric model at tree level. 
In this model, one of the two Higgs doublets couples to 
up- type quarks, while the other doublet couples to down- 
type quarks and leptons. 

The contributions of the charged Higgs to B — >• 
D^*'t~V t decays can be encapsulated in the scalar helic- 
ity amplitude in the following way [U, [13] : 



H 



2HDM 



H 



SM 



tan 2 /3 



m 2 H± 1 =F m c /m b 



(19) 



Here, tan/3 is the ratio of the vacuum expectation values 
of the two Higgs doublets, m#± is the mass of the charged 
Higgs, and m c /m b = 0.215 ±0.027 2l| is the ratio of the 
c- and 6-quark masses at a common mass scale. The 
negative sign in Eq. [19] applies to B — » Dt~v t decays 
and the positive sign applies to B — > D*t~V t decays. 
This expression is accurate to 1% for m H + larger than 
15 GeV. The region for m H + < 15 GeV has already been 
excluded by B — > X s ^/ measurements [22j|. 

The tan(3/m H + dependence of the ratios TZ(D^) in 
the type II 2HDM can be studied by substituting ff 2HDM 
for Hf M in Eq. [SJ Given that charged Higgs bosons are 
not expected to contribute significantly to B —> D^l~Vi 
decays, 7?.(D^*^)2hdm can be described by a parabola in 
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TABLE I. Dependence of TZ(D ( * ) ) on tan/3 /m H+ in the 
2HDM according to Eq.[20]for B -)• Dt~V t and B D*t~V t 
decays: the values of 1Z(D^), the parameters A and B with 
their uncertainties, and correlations C. 





B -> Dt~V t 


B -> D*t~V t 




0.297 ± 0.017 


0.252 ± 0.003 


A dW (GeV 2 ) 


-3.25 ± 0.32 


-0.230 ± 0.029 


B DM (GeV 4 ) 


16.9 ± 2.0 


0.643 ± 0.085 


c(h(d^) sm ,a dM ) 


-0.928 


-0.946 


C(TZ(D^) S m,B d( , } ) 


0.789 


0.904 


C{A Dlt) ,B DM ) 


-0.957 


-0.985 



the variable t&n 2 [3/m 2 H+ , 



^(P W )2HDM 



SM" 



tan 2 /3 



-B 



tan 4 /? 



2 i "uvi 4 



(20) 

Table U lists the values of A d m and B D („), which are 
determined by averaging over B° and B~ decays. The 
uncertainty estimation includes the uncertainties on the 
mass ratio m c /mf, and the FF parameters, as well as their 
correlations. 

Due to the destructive interference between the SM 
and 2HDM amplitudes in Eq. [TH1 charged Higgs con- 
tributions depress the ratios TZ(D^) for low values of 
tan/3/m#+. For larger values of tan/3/m#+, the Higgs 
contributions dominate and Ti-(D) and 1Z(D*) increase 
rapidly. As the coefficients of Table Q] show, the 2HDM 
impact is expected to be larger for 1Z(D) than for 1Z(D*). 
This is because charged Higgs contributions only affect 
the scalar amplitude P 2HDM , but B — > D*t~V t decays 
also receive contributions from H± , diluting the effect on 
the total rate. 

Figure [3] shows the impact of the 2HDM on the q 2 spec- 
trum. Given that the B and D mesons have spin J = 0, 
the SM decays B — > DW* —> Dtv proceed via P-wave 
for Jw = 1, and via S-w&ve for J\y* = 0. For the P- 
wave decay, which accounts for about 96% of the total 
amplitude, the decay rate receives an additional factor 
|pj,| 2 , which suppresses the q 2 spectrum at high values. 
Since charged Higgs bosons have Jh = 0, their contribu- 
tions proceed via 5-wave, and, thus, have a larger average 
q 2 than the SM contributions. As a result, for low values 
of tan/3/m#+ where the negative interference depresses 
i/ 2HDM , the q 2 spectrum shifts to lower values. For large 
values of taxtft /m H + , the Higgs contributions dominate 
the decay rate and the average q 2 significantly exceeds 
that of the SM. 

The situation is different for B — > D*t~v t decays be- 
cause the D* meson has spin Jd» = 1. The SM decays 
can proceed via S, P, or D-waves, while the decay via 
an intermediate Higgs boson must proceed via P-wave, 
suppressing the rate at high q 2 . 

When searching for charged Higgs contributions, it is 
important to account for the changes in the q 2 spectrum. 
This distribution has a significant impact on the analysis 
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FIG. 3. (Color online). Predicted q 2 distributions for (a) 
B —¥ Dt~T>t and (b) B — > D*t~V t decays for different values 
of tan/3/m H +. All curves are normalized to unit area. 



q 2 and m 2 ^, one of 



due to the close relation between 
the fit variables. 

Charged Higgs contributions also affect the \p* t \ dis- 
tribution. Given the spin of the Higgs boson and the 
positive hclicity (right-handedness) of the anti-neutrino, 
the decays H~ — > t~V t always produce t~ leptons with 
positive hclicities (A T = +). As a result, the fraction of 
right-handed leptons produced in B — >• D(*'t~V t de- 
cays changes from 30% in the SM, to close to 100% when 
the 2HDM contributions dominate. 

The lepton spectrum of polarized r — > l^viv T de- 
cays is well known (23[. For r~ leptons with A T - = — , 
the l~ is emitted preferentially in the r~ direction, while 
the opposite is true for positive hclicities. In the B rest 
frame, leptons of a certain momentum in the t~ rest 
frame have larger momentum if they are emitted in the 
direction of the r~ momentum than in the opposite di- 
rection. As a result, the \p\\ spectrum for SM decays is 
harder than for Higgs dominated decays. For low val- 
ues of tan/3 /m H + for which the destructive interference 
depresses the B — > D(*>t~V t rate, the proportion of left- 
handed t~ leptons increases, and therefore, the \p\ \ spec- 
trum is harder than in the SM. 



III. DATA SAMPLE, DETECTOR AND 
SIMULATION 

A. Data Sample 

This analysis is based on the full data sample recorded 
with the BAR4R detector [24[ at the PEP-II energy- 
asymmetric e + e~ storage rings [25j . It operated at a 
center-of-mass (cm.) energy of 10.58 GeV, equal to the 
mass of the T(4S) resonance. This resonance decays al- 
most exclusively to BB pairs. The collected data sample 
of 471 million T(4S) — > BB events (on-peak data), corre- 
sponds to an integrated luminosity of 426 fb _1 [26|]. To 
study continuum background, an additional sample of 
40fb _1 (off-peak data) was recorded approximately 40 
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McV below the T(AS) resonance, i.e., below the thresh- 
old for BB production. 



B. The BABAR Detector and Single Particle 
Reconstruction 

The BABAR detector and event reconstruction have 
been described in detail elsewhere J24[. The momentum 
and angles of charged particles were measured in a track- 
ing system consisting of a 5-laycr, double-sided silicon- 
strip detector (SVT) and a 40-layer, small-cell drift cham- 
ber (DCH) filled with a helium-isobutane gas mixture. 
Charged particles of different masses were distinguished 
by their ionization energy loss in the tracking devices and 
by a ring- imaging Cerenkov detector (DIRC). A finely 
segmented CsI(Tl) calorimeter (EMC) measured the en- 
ergy and position of electromagnetic showers generated 
by electrons and photons. The EMC was surrounded 
by a superconducting solenoid providing a 1.5-T mag- 
netic field and by a segmented flux return with a hexag- 
onal barrel section and two endcaps. The steel of the 
flux return was instrumented (IFR) with resistive plate 
chambers and limited streamer tubes to detect particles 
penetrating the magnet coil and steel. 

Within the polar angle acceptance of the SVT and 
DCH (0.4 < diab < 2.6) the efficiency for the recon- 
struction of charged particles exceeds 99% for momenta 
above 1 GeV. For low momentum pions, especially from 
D* + —> D°n + decays, the efficiency drops to about 90% 
at 0.4 GeV and to 50% at 0.1 GeV. 

The electron and muon identification efficiencies and 
the probabilities to misidentify a pion, a kaon, or a proton 
as an electron or a muon are measured as a function of 
the laboratory momentum and angles using high-purity 
data samples. 

Electrons are separated from charged hadrons primar- 
ily on the basis of the ratio of the energy deposited in the 
EMC to the track momentum. A special algorithm has 
been developed to identify photons from bremsstrahlung 
in the inner detector, and to correct the electron mo- 
mentum for the energy loss. Within the polar angle ac- 
ceptance, the average electron efficiency for laboratory 
momenta above 0.5 GeV is 97%, largely independent of 
momentum. The average pion misidentification rate is 
less than 0.5%. 

Muon identification relies on a new multivariate algo- 
rithm that significantly increases the reconstruction effi- 
ciency at low muon momenta, \p \ < 1 GeV. This algo- 
rithm combines information on the measured DCH track, 
the track segments in the IFR, and the energy deposited 
in the EMC. The average muon efficiency is close to 90% 
independent of momentum, except in the forward end- 
cap, where it decreases for laboratory momenta below 
1 GeV. The average pion misidentification rate is about 
2% above 1.2 GeV, rising at lower momenta and reaching 
a maximum of 9% at 0.8 GeV. 

By choosing a fairly loose selection of charged leptons 



and taking advantage of improved PID algorithms, we 
increased the lepton efficiencies by 6% for electrons and 
50% for muons compared to the previous BABAR analy- 
sis [3. 

Charged kaons are identified up to 4 GeV on the basis 
of information from the DIRC, SVT, and DCH. The ef- 
ficiency exceeds 80% over most of the momentum range 
and varies with polar angle. The probability that a pion 
is misidentified as a kaon is close to 2%, varying by about 
1% as a function of momentum and polar angle. 

The decays K® — > 7r + 7r~ are reconstructed as pairs 
of tracks of opposite charge originating from a displaced 
vertex. The invariant mass of the pair m T7r is required to 
be in the range G [0.491, 0.506] GeV. No attempt is 
made to identify interactions of in the EMC or IFR. 

To remove beam-generated background in the EMC 
and electronic noise, photon candidates arc required 
to have a minimum energy of 30 MeV and a shower 
shape that is consistent with that of an electromagnetic 
shower. Neutral pions are reconstructed from pairs of 
photon candidates with an invariant mass in the range 
m 77 e [120, 150] MeV. 



C. Monte Carlo Simulation 

1. Simulated Samples 

This analysis relics on Monte Carlo (MC) techniques to 
simulate the production and decay of continuum and BB 
events. The simulation is based on the EvtGen generator 
[27| . The qq fragmentation is performed by Jetset f28j . 
and the detector response by Geant4 [2^]. Radiative ef- 
fects such as bremsstrahlung in the detector material and 
initial-state and final-state radiation [3(| are included. 

We derive predictions for the distributions and efficien- 
cies of the signal and backgrounds from the simulation. 
The size of the simulated sample of generic BB events 
exceeds that of the BB data sample by about a factor of 
ten, while the sample for qq events corresponds to twice 
the size of the off-peak data sample. We assume that the 
T(AS) resonance decays exclusively to BB pairs and use 
recent measurements of branching fractions fl2j for all 
produced particles. The impact of their uncertainties on 
the final results is assessed as a systematic uncertainty. 

Information extracted from studies of selected data 
control samples is used to improve the accuracy of the 
simulation. Specifically, we reweight simulated events to 
account for small differences observed in comparisons of 
data and simulation fSec. fV|) . 



2. Implementation of the Form Factor Parameterizations 

For reasons of simplicity, the simulation of B — > D£~Vg 
and B -> D [ -*^t^V t decays is based on the ISGW2 
model [3l|, and B —> D*£~V( decays are generated using 
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an HQET-bascd parameterization [32j . A change to a dif- 
ferent FF parameterization is implemented by reweight- 
ing the generated events with the weights 



/ A4(g 2 ,0»)HQET 

^ M{q\e t ) MC 



B 



MC 



HQET 



(21) 



Here, A4(q 2 , 6\)hqet refers to the matrix element for 
the FF parameterizations described in Sees. Ill A II and 
III A 21 and M(q 2 , $i)MC is the matrix element employed 
in the MC generation. The matrix element of decays 
involving the scalar D meson depends on one angular 
variable, the lepton helicity angle 0£, with I = e, /i, r. In 
addition to 9e, the matrix element of decays involving 
the vector meson D* is sensitive to two additional angu- 
lar variables describing the D* decay. The ratio of the 
branching fractions Bmc/Bhqet ensures that the sum of 
all weights equals the number of generated events. 

In the SM, this reweighting results in a small shift of 
the q 2 distribution to higher values, while the changes in 
the helicity angle 9 T and the r polarization are negligi- 
ble. Therefore, the distributions of the secondary charged 
lepton are not affected. 

In the presence of a charged Higgs boson, however, 
the t polarization can change substantially, affecting the 
momentum of the secondary lepton £ originating from the 
r — > lv(i> T decays. We account for the potential presence 
of a charged Higgs of 2HDM type II by reweighting the 
simulation with the following weights, 



w 2 HDM(q 2 ,9i, \Pt\ 



M{q 2 A) 



2HDM 



X(q 2 ,^)MC 



r(|p||) 



2HDM 



B 



MC 



MC 



B 



2HDM 



(22) 



where Qi refers again to the angular variables. The second 
factor represents the ratio of the \p\ | distributions r(|p||) 
in the 2HDM parameterization and in the MC simula- 
tion. This factorization is necessary because in the MC 
generation the polarization is handled in a probabilistic 
manner, so it cannot be corrected on an event- per- event 
basis. It is only applicable if \p\ | is uncorrelated with 
q 2 and the angular variables, which is largely the case. 
In some regions of phase space, the 2HDM weights have 
a much larger dispersion than the weights applied in the 
SM reweighting, leading to larger statistical uncertainties 
for the simulation of the Higgs boson contributions. 



arc not included in the nominal fit for lack of reliable 
measurements. 

To estimate the rate of B D**tv t decays, we rely 
on ratios of the available phase space $, 



Tl{D**) 



B(B -> D**t- 



$(B -> D**t- 



B(B -> D**£-V e ) <f>(B 



D**£~vA 



(23) 

The value of this ratio depends on the mass of the D** 
state involved in the B — > D**(t~ /£~)v decay. We use 
the largest of the four possible choices, 1Z(D**) = 0.18. 
Possible contributions from non-rcsonant B — > 
decays and semileptonic decays involving 
higher-mass excited charm mesons are not included in 
the nominal fit, and will be treated as a systematic un- 
certainty. 



IV. EVENT SELECTION 

The event selection proceeds in two steps. First, we 
select BB events in which one of the B mesons, the 
Stag, is fully reconstructed in a hadronic decay, while the 
other B meson decays semileptonically. To increase the 
event selection efficiency compared to earlier analyses, 
we have added more decay chains to the -Btag selection 
and have chosen a looser charged lepton selection. This 
leads to significantly higher backgrounds, primarily com- 
binatorial background from BB and continuum events, 
and charge-crossfeed events. Charge-crossfeed events are 
B — > (t~ I '£~)v decays in which the charge of the re- 
constructed Btag and mesons are wrong, primarily 
because of an incorrectly assigned low-momentum 7T . 

Semileptonic decays to higher mass charm mesons have 
a signature similar to that of signal events and their com- 
position is not well measured. This background is fitted 
in selected control samples that are enriched with these- 
decays. 

As the second step in the event selection, we introduce 
kinematic criteria that increase the fraction of selected 
signal events with respect to normalization and back- 
ground decays. We also apply a multivariate algorithm 
to further improve the signal-to-background ratio. 



Selection of Events with a Stag and a 
Semileptonic B Decay 



3. Simulation of B — > D** (r jl )v decays 

By D** we refer to excited charm resonances heavier 
than the D* meson. We include in the simulation the 
B -> D**t~V t and B — > D**£~Vi decays that involve 
the four D** states with L = 1 that have been measured 
. This simulation takes into account their helicitics 
33j and the following decay modes: Z?q,£>2 — > Dtt and 
D[,Di,D* -> D*tt. Three-body decays D** -> D^tttt 



T{AS) — > BB events are tagged by the hadronic decay 
of one of the B mesons. We use a semi-exclusive algo- 
rithm which includes additional -Btag decay chains and 
enhances the efficiency by a factor of 2 compared to the 
earlier version employed by BABAR |l4[. We look for de- 
cays of the type Stag — > SX^, where S refers to a seed 
meson and is a charged state comprising of up to five 
hadrons, pions or kaons, among them up to two neutral 
mesons, 7r° or The seed mesons, D, D* , D s . D*, and 
J/0, are reconstructed in 56 decay modes. As a result, 
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the -Btag is reconstructed in 1,680 different decay chains, 
which arc further subdivided into 2,968 kinematic modes. 

To isolate the true tag decays from combinatorial back- 
ground, we use two kinematic variables: the energy sub- 
stituted mass rriES = \J ^ cam — Pt ag and the energy dif- 
ference AE = Stag - -Bboam- Here p tag and E t&g refer to 
the cm. momentum and energy of the -B tag , and -Ebcam is 
the cm. energy of a single beam particle. These variables 
make optimum use of the precisely known energies of the 
colliding beams. For correctly reconstructed B decays, 
the m,ES distribution is centered at the B- meson mass 
with a resolution of 2.5 MeV, while AE is centered at 
zero with a resolution of 18 MeV which is dominated by 
the detector resolution. We require dies > 5.27 GeV and 
\AE\ < 0.072 GeV. 

For each B tag candidate in a selected event, we look 
for the signature of the semileptonic decay of the sec- 
ond B meson, a D or D* meson and a charged lepton 
I. We combine charged Stag candidates with D^*'°£~ 
systems and neutral -B tag candidates with both D^ + £" 
and t + systems, where the inclusion of both charge 

combinations allows for neutral B mixing. We require 
all charged particles to be associated with the B tag D^£ 
candidate, but we allow for any number of additional 
photons in the event. 

The laboratory momentum of the electron or muon 
is required to exceed 300 MeV or 200 MeV, respec- 
tively. For D mesons, we reconstruct the following decay 
modes: D° K~tt+, K~K+, K'-k+tt , K-tt+tt~tt+ } 
K^tt+tt-, and D+ K~tt+tt+, K-tt+tt+tt , K°tt+ , 

K^TT+TT+TT-, K° S TT + TT°, K°K+ , with K° TT+TT-. The 

reconstructed invariant mass of D candidates is required 
to be consistent with the nominal D mass to within four 
standard deviations (a). The combined reconstructed 
branching fractions are 35.8% and 27.3% for D° and 
D + , respectively. We identify D* mesons by their decays 
D*+ D q tt+,D+tt q } and D*° D Q ir a ,D a j. For these 
decays, the cm. momentum of the pion or the cm. en- 
ergy of the photon are required to be less than 400 MeV. 
Furthermore, the mass difference Am = m(D*) ~ m(D) 
is required to differ by less than 4<r from the expected 
value fl2| . 

To further reduce the combinatorial background, we 
perform a kinematic fit to the event, constraining tracks 
of secondary charged particles to the appropriate B, 
D^*\ or K® decay vertices. The fit also constrains the 
reconstructed masses of the D, D*, and mesons to 
their nominal values. The vertex of the T(AS) —> BB 
decay has to be compatible with a beam-beam interac- 
tion. Candidates for which this fit does not converge are 
rejected. The m^ iss resolution improves by about 25% 
and becomes more symmetric for the remaining candi- 
dates. 

To select a single BB candidate, we determine -Eextra = 
J2iE] , the sum of the energies of all photons that are not 
associated with the reconstructed BB pair. We only in- 
clude photons of more than 50 MeV, thereby eliminating 



about 99% of the beam-generated background. We re- 
tain the candidate with the lowest value of -E oxtra , and if 
more than one candidate survives, we select the one with 
the smallest \AE\. This procedure preferentially selects 
D*£ candidates over D£ candidates. Thus, we reduce the 
fraction of misreconstructed events with a D* — > D^/j) 
decay for which the pion or photon is not properly as- 
signed to the D* meson. 

As a consequence of the rather loose lepton selection 
criteria and the addition of decay modes with multiple 
neutral pions and K® for the Btag selection, the number 
of BtzgD^C candidates per event is very large. To ad- 
dress this problem, we identify the -B tag decay modes that 
contribute primarily to the combinatorial background. 
Specifically, we determine for each of the 2,968 kine- 
matic modes i?tc, the fraction of events for which all 
charged particles in the B tag final state are correctly re- 
constructed and associated with the tag decay. This as- 
sessment is based on a large sample of simulated BB 
events equivalent to 700 fb -1 . We observe that for de- 
cay chains with low multiplicity final states and no neu- 
tral hadrons the signal-to-background ratio {S/B) is very 
high. For instance, for the -Btag ~~ * ^AH - * A t+ / i_ )^' _ 
decay, we obtain S/B = 316/79, whereas for the de- 
cay J5{* ag — > £)"(—>• A'°7r~ )7T + 7r + 7r + 7r~7r~ this ratio is 
S/B = 20/145. For this decay mode, typically 3.5 of 
the 8 i?t ag final state particles are incorrectly associ- 
ated with the second B decay in the event or otherwise 
misidentified. Based on this study, we only retain -B tag 
decay chains with R tc > 0.3. With this criterion, we re- 
move 2100 -Bt a g kinematic modes, eliminate 2/3 of the 
combinatorial background, and retain 85% of the signal 
B — > D^*'t~V t decays. Thanks to this procedure, the av- 
erage number of candidates per event before single can- 
didate selection is reduced to 1.8 for the D°£ and D + £ 
samples, and 3.1 and 4.8 for the D*°£ and D* + £ samples, 
respectively. 

B. Selection of the D^-K°t Control Samples 

To constrain the B — > D**(t~ / t~)v background, we 
select four D^tt°£ control samples, identical to the D^*>£ 
samples except for an additional reconstructed ir°. The 
7T° is selected in the mass range m 77 G [120, 150] MeV. 
Decays of the form B — s- D^it£v peak at m^ liss = in 
these samples. As a result, we can extract their yields to- 
gether with the signal and normalization yields by fitting 
the D^i and D^tt°£ samples simultaneously. 

More than half of the events in these control samples 
originate from continuum e + e~ — > qq{j) events. Since 
the fragmentation of light quarks leads to a two-jet event 
topology, this background is very effectively suppressed 
by the requirement | cos A6* t hrust | < 0.8, where A# t hrust 
is the angle between the thrust axes of the B tag and of 
the rest of the event. Since B mesons originating from 
T(AS) decays are produced just above threshold, their 
final state particles are emitted almost isotropically, and, 
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decays. The event sample with q 2 < 4 GeV is domi- 
nated by B — > D^£~T>i and serves as a very clean data 
sample for comparisons with the MC simulation. To re- 
ject background from hadronic B decays in which a pion 
is misidentified as muon, we require |p m j SS | > 200 MeV, 
where |p m j SS | is the missing momentum in the cm. frame. 

To further improve the separation of well- 
reconstructed signal and normalization decays from 
various backgrounds, we employ a boosted decision 
tree (BDT) multivariate method [IT 
relics on simple classifiers which 
and background regions by using 



| . This method 
determine signal 
binary selections 



on various input distributions. For each of the four 



select signal and 
-> D**(t-/1~)V 



D^*H samples, we train a BDT to 
normalization events and reject B 
and charge cross- feed backgrounds. Each BDT selector 
relics on the simulated distributions of the following 
variables: (a) -E GX tra; (b) AE; (c) the reconstructed mass 
of the signal D^*> meson; (d) the mass difference for the 
reconstructed signal D*: Am = m(Dir) — m(D); (e) 
the reconstructed mass of the seed meson of the Stag 5 
(f) the mass difference for a D* originating from the 
Stag, Am tag = m(AagTr) - m(D tag ); (g) the charged 
particle multiplicity of the Stag candidate; and (h) 
cos A6>thrust • The input distributions for one of the BDT 
selectors are shown in Fig. [4] For the D^n°£ samples, 
we use similar BDT selectors that are trained to reject 
continuum, £)(*) (£ / 't)v , and other BB background. 
After the BDT requirements are applied, the fraction 
of events attributed to signal in the m^ iss > 1.5 GeV 2 
region, which excludes most of the normalization decays, 
increases from 2% to 39%. The background remaining 
in that region is composed of normalization events 
(10%), continuum (19%), D**{£/r)v events (13%), and 
other BB events (19%), primarily from B -> D^D^*^ 
decays with Df — > t + i> t . 



FIG. 4. (Color online). Input variables for the BDT selector 
trained on the D*°£ sample. Histograms are normalized to 
1000 entries. 

therefore, the cosA0 t hrust distribution is uniform. As 
a result, the loss of B —> D**(t~ / £~)v decays due to 
this restriction is significantly smaller than the amount 
of continuum events rejected. 

C. Optimization of the Signal Selection 

We introduce criteria that discriminate signal from 
background, and also differentiate between signal B — > 
D(**>t~V t and B — > D^£~T>i decays. For semileptonic 
decays the minimum momentum transfer is largely de- 
termined by the mass of the charged lepton. For decays 
involving r leptons, = m 2 . ~ 3.16 GeV 2 . Thus the 
selection q 2 > 4 GeV 2 retains 98% of the B_-+ D^t~V t 
decays and rejects more that 30% of the B — s- D^*>t~Vi 



V. CORRECTION AND VALIDATION OF THE 
MC SIMULATION 

The simulation of the full reconstruction of high- 
multiplicity events, including the veto of events with ex- 
tra tracks or higher values of -E e xtra is a rather challenging 
task. To validate the simulation, we compare simulated 
distributions with data control samples, and, when neces- 
sary, correct the MC simulations for the observed differ- 
ences. The figures shown in this section combine events 
from all four channels (D°l, D*°£, D+£, and D*+£); the 
observed differences are similar in the individual samples. 

The control samples are selected to have little or no 
contamination from signal decays. Specifically we select, 

• Continuum events: off-peak data. 

• Normalization decays: q 2 < 4 GeV 2 . 

• Combinatorial BB and continuum backgrounds: 
5.20 < m ES < 5.26 GeV 
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FIG. 5. (Color online). Comparison of data control samples (data points) with MC simulated samples (histograms) of the |pj| 
distributions for (a) off-peak data prior to |p|| reweighting, (b) the intermediate -Eextra sample prior to |p|| reweighting, and 
(c) the intermediate .Eextra sample after |p^| reweighting; (d) the -Eextra distribution for the combinatorial background; and the 
771ES distributions for (e) the intermediate i? ex tra sample, and (f) the low -E cx tra sample. The results are shown for the four 
D^i samples combined. 



• Incorrectly reconstructed events: events in three 
Sextra intervals, high (1.2 < E extla < 2.4 GeV), 
intermediate (0.5 < -E C xtra < 1.2 GeV), and low 
(£ oxtra < 0.5 GeV for events that fail the BDT 
selection). N.B. the BDT selection results in the 
elimination of all events with i? oxtra > 0.4 GeV. 

The off-peak data sample confirms the to^^ distribu- 
tion of simulated continuum events, but shows discrepan- 
cies in the \p* t | spectrum and overall normalization of the 
simulation [Fig. Ufa)]. These features are also observed 
in other control samples, such as on-pcak data with high 
£-oxtra [Fig.E^b)]. We correct the simulated |p^| spectrum 
and yield of the continuum contribution by reweighting it 
to match off-peak data, on an event-by-event basis. Af- 
ter this correction, the \p* t \ distributions of the expected 
backgrounds agree well in independent control samples 
down to low lepton momenta where the misidentification 
rates are significant [Fig. Etc)]. Wc observe that in the 
high -E ex tra region, the simulation exceeds data yield by 
(1.3±0.5)%. This small excess is corrected by decreasing 
the expected BB background yield by (4.3 ± 1.9)%. Af- 
ter this correction, the simulation provides accurate yield 
predictions for the backgrounds at intermediate and high 
-S'extra- For instance, the ratio of the expected to observed 
yield of events with m^ lss > 1.5 GeV 2 is 0.998 ± 0.006. 
The m. 2 ^^ distributions of the continuum and BB back- 
grounds are described well in all control samples. 

The region of low -Eextraj which includes the signal re- 
gion, is more difficult to model, primarily due to low 



energy photons and K® mesons interacting in the EMC. 
Figure Eld) shows that the data in the ragg sideband 
agree well with the combinatorial background predictions 
for -Ecxtra > 0.5 GeV, but are underestimated for low 
^cxtra- This, and small differences in the other BDT in- 
put distributions, result in a underestimation of the com- 
binatorial background when the BDT requirements are 
applied. Based on the 5.20 < toes < 5.26 GeV sideband, 
we find scale factors of 1.099±0.019 and 1.047±0.034 for 
the combinatorial background in the D£ and D*£ sam- 
ples, respectively. The uncertainties are given by the 
statistics of the data and simulated samples. The ra- 
tio between the observed and expected Toes distribu- 
tion is independent of -E C xtra [Figs. EJe,f)], so we ap- 
ply these corrections to the continuum and BB back- 
grounds in the signal region. The same correction is 
applied to B —> D**(t~ / l~)v events, which cannot be 
easily isolated, because their simulated -Eextra distribu- 
tions are very similar to those of combinatorial back- 
ground. These corrections affect the fixed BB and con- 
tinuum yields in the fit, as well as the relative efficiency 
of B -> D**(t-/£-)V events in the D^t and D^TT°e 
samples. As a result, these corrections are the source of 
the dominant systematic uncertainties. 

Relying on the q 2 < 4 GeV 2 control sample, where 
B — > D(*H~T>t decays account for 96% of the events, 
we correct the E ext r& distribution and an 8.5% overesti- 
mation of the simulated normalization events. We apply 
the same correction to simulated signal events which are 
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expected to have a similar -E C xtra distribution. This pro- 
cedure does not affect the relative efficiency of signal to 
normalization events, so it has a very small impact on 
the TZ(D^) measurements . 

We use the same q 2 < 4 GeV 2 control sample to com- 
pare and validate the \p* t \ distributions of B — >• D^£~Ve 
events. We observe that the TO^ iss resolution of the nar- 
row peaks at mjju S8 = is slightly underestimated by 
the simulation. This effect is corrected by convolving the 
simulated distributions with a Gaussian resolution func- 
tion, for which the width is adjusted by iteration. 



VI. FIT PROCEDURE AND RESULTS 
A. Overview 

We extract the signal and normalization yields from 
an extended, unbinned maximum-likelihood fit to two- 
dimensional wjuss - distributions. The fit is per- 
formed simultaneously to the four D^*'£ samples and the 
four D^tt°£ samples. The distribution of each D^£ and 
D^7r°£ sample is fit to the sum of eight or six contribu- 
tions, respectively. Each of the 4x8 + 4x6 = 56 con- 
tributions is described by a probability density function 
(PDF). Their relative scale factor determines the number 
of events from each source. Tables HTl and ITO1 summarize 
the contributions to the fit for the four D^£ sample and 
the four D^*>tt°£ samples. These tables also list the rel- 
ative yield for each contribution as estimated from MC 
simulation (for SM signal), and specify whether the yield 
is free, fixed, or constrained in the fit. 

We introduce the following notation to uniquely iden- 
tify each contribution to the fit: source => sample. For 
instance, D*°tv =>■ D*°£ refers to signal D*®tv de- 
cays that are correctly reconstructed in the D*°£ sam- 
ple, while D*°tv =>■ D°£ refers to the same decays, but 
incorrectly reconstructed in the D°£ sample. We refer 
to the latter as feed-down. Contributions of the form 
D{t/£)v =>> D*(t/£) and D^{t/£)v =$> D**{t/£) are re- 
ferred to as feed-up. 

The contributions from the continuum, BB, and cross- 
feed backgrounds, with the exception of BB background 
in the D^tt°£ samples, are fixed to the yields determined 
by MC simulation after small adjustments based on data 
control regions. The yields of the remaining 36 contri- 
butions are determined in the fit. Some of these con- 
tributions share the same source and therefore the ra- 
tio of their yields is constrained to the expected value, 
e.g., D* a rv =S> D*°£ and D* a rv =S> D°£. Of special im- 
portance are the constraints linking the D**(£/t)v yields 
in the D^*>i samples [N D ,,^ D (,)) to the yields in the 
D^n°£ samples (N D «,^ DM „ ), 
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D • 



(24) 



is equivalent to the ratio of the D**(£/t)v reconstruction 
efficiencies for the two samples. 

Taking into account the constraints imposed on event 
yields from a common source, there are 22 free parame- 
ters in the standard fit, as listed in Table ITVl In addition, 
we perform a fit in which we impose the isospin relations 
K{D°) = K(D+) = 11(D) and U(D*°) = U(D*+) = 
1Z(D*). We choose not to impose isospin relations for 
the D(*>tt°£ samples. Consequently, this fit has a total 
of 17 free parameters. 

The following inputs are updated by iterating the fit: 

• The eight D^(£/t)u => D^n £ PDFs are recal- 
culated taking into account the fitted D^£v and 

tv contributions to the D^£ samples. 

• The fixed charge cross-feed yields are updated 
based on the deviation of the fitted D^£v yields 
from the expected values. 

• The continuum, BB, and D**(£/t)v background 
corrections are recalculated. They have a slight de- 
pendence on the fitted D^£v events because some 
of these events extend into the toes sideband. 



The correction to the 



resolution of the nor- 



malization contributions is readjusted. 

• The two feed-down constraints for D*tv are up- 
dated using the fitted feed-down constraints for the 
normalization contributions in the following way: 
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D> 



■D*e 
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MC 



N D . 
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(25) 



MC 



Given that these constraints share the same source, /d** 



The iterations continue until the change on the values of 
TZ(D^) is less than 0.01%. The update of the feed-down 
rates has a significant impact on the fits to the D° and 
D + samples because of the large signal feed-down. The 
other iterative updates have only a marginal impact. 



B. Probability Density Functions and Validation 

The fit relies on 56 PDFs, which are derived from MC 
samples of continuum and BB events equivalent to 2 and 
9 times the size of the data sample, respectively. The 
two-dimensional ?Ti^ iss -|p|| distributions for each of the 
56 contributions to the fit are described by smooth non- 
parametric kernel estimators (35|. These estimators en- 
ter a two-dimensional Gaussian function centered at the 
m miss an d \p*i \ values of each simulated event. The width 
of the Gaussian function determines the smoothness of 
the PDF. We find the optimum level of global smooth- 
ing with a cross-validation algorithm J36|. For PDFs 
that have variations in shape that require more than one 
level of smoothing, we combine estimators with different 
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TABLE II. Contributions to the four D^*'£ samples. The expected relative abundance of events in each data sample is 
represented by / exp . The columns labeled Yield indicate whether the contribution is free in the fit, fixed, or linked to another 
component through a cross-feed constraint. The charged cross-feed components, marked with Fix./It., are fixed in the fit, but 
updated in the iterative process. 





D° 


£ 


D* 


°e 


D+ 


I 


D*~ 


v £ 


Source 


/cxp (%) 


Yield 


/cxp (%) 


Yield 


/cxp (%) 


Yield 


/cxp (%) 


Yield 


D^rv signal 


2.6 


Free 


4.9 


Free 


4.3 


Free 


5.0 


Free 


D^*'tv signal feed-down/up 


2.8 


D*°£ 


0.4 


dH 


1.8 


D* + £ 


0.1 


D+£ 


D^£v normalization 


24.5 


Free 


80.7 


Free 


37.3 


Free 


88.0 


Free 


D^lv norm, feed-down/up 


53.5 


Free 


2.7 


D°£ 


35.0 


Free 


0.3 


D+£ 


D** (£/t)u background 


4.3 


D°tv°£ 


3.6 


D*°ir°£ 


6.6 


D+t: £ 


3.0 


D*+-k°1 


Cross-feed background 


3.8 


Fix./It. 


1.3 


Fix./It. 


2.1 


Fix./It. 


0.4 


Fix./It. 


BB background 


4.1 


Fixed 


3.7 


Fixed 


7.1 


Fixed 


2.8 


Fixed 


Continuum background 


4.4 


Fixed 


2.6 


Fixed 


5.9 


Fixed 


0.5 


Fixed 



TABLE III. Contributions to the four D^*'it £ samples. The expected relative abundance of events in each data sample is 
represented by / exp . The columns labeled Yield indicate whether the contribution is free in the fit, fixed, or linked to another 
component through a cross-feed constraint. The D(£/t)u components are linked to the D*{£/r)u components, and the cross- 
feed constraint is updated in the iteration. The charged cross- feed components, marked with Fix./It., axe fixed in the fit, but 
updated in the iterative process. 



Source 


D°TT°£ 

/cxp (%) Yield 


D*°ir°£ 
/cxp (%) Yield 


D+n°£ 
/cxp (%) Yield 


D*+tv°£ 
/cxp (%) Yield 


D**(£/t)v background 
D*(£/r)u feed-up 
D(£/t)v feed-up 
Cross-feed background 
BB background 
Continuum background 


20.1 Free 
19.1 Free 
6.4 D°tv°£ 
4.9 Fix./It. 
28.4 Free 
21.0 Fixed 


16.4 Free 
20.6 Free 
2.3 D*°ir°£ 
3.6 Fix./It. 
36.4 Free 
20.8 Fixed 


19.9 Free 
10.0 Free 
4.7 D + tv°£ 
4.4 Fix./It. 
38.7 Free 
22.2 Fixed 


22.1 Free 

25.2 Free 
0.8 D* + tt°£ 
2.5 Fix./It. 

37.4 Free 
12.0 Fixed 



TABLE IV. Number of free parameters in the isospin- 
unconstrained (N un ) and constrained (N cons ) fits. 



Sample Contribution 7V un JV t 





£ 


D^*'rv signal 


4 


2 


D (* 


£ 


D^*'£v normalization 


4 


2 


D (* 


£ 


D*£v norm, feed-down 


2 


1 


D (* 


7T°£ 


D"(£/t)v background 


4 


4 


D (* 


7Y°£ 


D^*'£v norm, feed-up 


4 


4 


D (* 


7T°£ 


BB background 


4 


4 



Gaussian widths in up to four areas in the n^iss - !? 5 !! 
space. For instance, we use different levels of smoothing 
in the D*®lv =>■ D*°£ contribution for the narrow peak 
at m^ liss = and the smooth m^ iss tail that extends 
up to 7GeV 2 . Figure [5] shows one-dimensional projec- 
tions of five two-dimensional PDFs. The bands indicate 
the statistical uncertainty on the PDFs estimated with a 
bootstrap algorithm (36j . 

The m,^ liss distributions of signal and normalization are 
very distinct due to the different number of neutrinos 



in the final state. The m^ liss distributions of the back- 
grounds resemble those of the signal, and therefore these 
contributions to the fit are either fixed or constrained by 
the D^n°£ samples. 



To validate the PDFs and the fit procedure, we divide 
the large sample of simulated BB events into two: sam- 
ple A with about 3.3 x 10 9 BB events, and sample B with 
9.4 x 10 s BB events. We determine the PDFs with sam- 
ple A, and create histograms by integrating the PDFs 
in bins of their ™^ iss and projections. We com- 
pare the resulting histograms with the events in sample 
A, and derive a \ 2 based on the statistical significance 
of the difference for each bin. The distribution of the 
corresponding p values for these PDFs is uniform, as ex- 
pected for an unbiased estimation. As another test, we 
extract the signal and normalization yields from fits to 
the events of sample B, using the PDFs obtained from 
sample A. Again, the results are compatible with an un- 
biased fit. Furthermore, we validate the fit procedure 
based on a large number of pseudo experiments gener- 
ated from these PDFs. Fits to these samples also show 
no bias in the extracted signal and normalization yields. 



17 




mLss (GeV 2 ) |pl| (GeV) 

FIG. 6. (Color online) Projections of the simulated m^iss 
and \p* t | distributions and the PDFs for the following contri- 
butions to the D°£ sample: (a), (b) D°tv; (c),_(d) D°lv; (e), 
(f) D*°£u; (g), (h) D"(£/t)u, and (i), (j) BB background. 
The light and dark blue (gray) bands mark the la and 2a en- 
velopes of the variations of the PDF projections due to their 
statistical uncertainty. 



C. Fit Results 

Figures [7] and [5] show the m^ iss and \p* t \ projections 
of the fits to the D^H samples. In Fig. the \p}\ 
projections do not include events with TO 2 liss > lGeV 2 , 
i.e., most of the signal events. In Fig. the vertical 
scale is enlarged and the horizontal axis is extended for 
the m^jgj, projection to reveal the signal and background 
contributions. The projections emphasize the signal 
events by excluding events with rn^ iss < 1 GeV 2 . Both 
figures demonstrate that the fit describes the data well 
and the observed differences are consistent with the sta- 
tistical and systematic uncertainties on the PDFs and 
the background contributions. 

Figure [5] shows the and projections of the fit 

to the four D^tt°£ samples. The narrow m 2 ^ peak is 
described well by the fit. It tightly constrains contribu- 
tions from B — > D^'niv decays, including the nonreso- 
nant D^tt states as well as decays of D** states, narrow 



or wide. There appears to be a small excess of events 
in the data for 1 < m 2 liss < 2 GeV 2 . This might be 
an indication for an underestimation of the D**(£/t)u 
background. The impact of this effect is assessed as a 
systematic uncertainty. 

The fit determines, for each signal decay mode, the 
number of signal events in the data sample, -/V s i g , and the 
corresponding number of normalization events, N nmm - 
We derive the ratios of branching fractions as 

ft(D(*>) = ^ sis £norm , (26) 

where £ s i g /£norm is the ratio of efficiencies (including 
the t branching fractions) taken from MC simula- 
tion. These relative efficiencies are larger for TZ-(D) than 
for 1Z(D*), because the q 2 > 4 GeV 2 requirement re- 
jects a larger fraction of B — > Dl~Vi decays than of 
B — > D*l~Vi decays, while keeping almost 100% of 
B — i> D^*'t~V t decays.. 

The results of the fits in terms of the number of events, 
the efficiency ratios, and TZ(D^) are listed in Table IVlnl 
for both the standard and the isospin-constraincd fits. 
Due to the large signal feed-down, there are significant 
negative correlations between the fits to the D£ and D*l 
samples. The statistical correlations are —0.59 for 7Z(D°) 
and K(D*°), -0.23 for K(D + ) and K(D*+), and -0.45 
for 11(D) and K(D*). 

VII. SYSTEMATIC UNCERTAINTIES 

Table [V] lists the systematic uncertainties considered 
in this analysis, as well as their correlations in the mea- 
surements of TZ(D) and 1Z(D*). We distinguish two kinds 
of uncertainties that affect the measurement of TZ(D^): 
additive uncertainties which impact the signal and back- 
ground yields and thereby the significance of the results, 
and multiplicative uncertainties that affect the £ s ig/£norm 
ratios and, thus, do not change the significance. The lim- 
ited size of the simulated signal and background samples 
impact both additive and multiplicative uncertainties. 

A. Additive uncertainties: 

Additive uncertainties affect the results of the fit. To 
asses their impact, we vary the source of uncertainty 1000 
times following a given distribution, and repeat the fit for 
each variation. We adopt as the uncertainty the standard 
deviation of the distribution of the resulting R(D^) val- 
ues. From this ensemble of fits, we also estimate the cor- 
relation between the uncertainties of 71(D) and 1Z(D*). 

1. PDF Estimation 

MC statistics: We employ a bootstrap algorithm [36[ 
to estimate the uncertainty due to the limited size of 
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FIG. 7. (Color online). Comparison of the m^ iss and |p|| distributions of the D^*'l samples (data points) with the projections 
of the results of the isospin-unconstrained fit (stacked colored distributions). The \p* t \ distributions show the normalization- 
enriched region with m^iss < 1 GeV 2 , thus excluding most of the signal events in these samples. 



the simulated event samples on which we base the 56 
PDFs. We generate 1000 samples of simulated events 
by sampling the original MC sample with replacement 
[37j . The PDFs are recalculated with each bootstrapped 
sample, and the fit is repeated for each set of PDFs. 
Figure [6] shows the la and 2a bands for the projections 
of five selected PDFs. The impact on the final result is 
4.4% for K(D) and 2.0% for K(D*). 

Form factors for B — s- D^*> (r~ / i~)V : We estimate 
the impact on the signal and normalization PDFs due 
to the uncertainties on the FF parameters, p 2 D , A, p 2 D , , 
i?o(l), -Ri(l), and i?2(l) ; taking into account their uncer- 
tainties and correlations. We recalculate the D^*'tv and 
D^iv PDFs with each set of 1000 Gaussian variations 
of the parameter values, and repeat the fit with each set 
of PDFs to determine the impact on TZ(D^). 

D** -> D ( *)(7r / 7 r ± ) fraction: The simulation of 
D**{1/t)v decays only includes the two-body decays 
D** -> D^n of the four L = 1 charm meson states. 
The ratio of D** -> D^tt decays to D** -> D^% ± 
decays which is fixed by isospin relations has a signifi- 



cant impact on the PDFs, because D** — > D^ir de- 
cays result in a sharply peaked m^ iss distribution for the 
D^7r°£ samples. The measured uncertainty on the tt° 
detection efficiency is 3%. We assume a 4% uncertainty 
to the probability that a low momentum charged pion 
from D** — > D^tt^ decays is misassigned to the B tag 
decay. Combining these two uncertainties, we arrive at 
an uncertainty on the relative proportion of the two-body 
decays of D** of 5%. We repeat the fit increasing and de- 
creasing this ratio by 5%, and adopt the largest variation 
of the isospin-constrained fit results as the systematic un- 
certainty. 

B — > D**t~T>i branching fractions: Since decays 
to the four D** states are combined in the B —> 
D(*\t~/£~)V samples, the PDFs depend on the rel- 
ative B — > D**i~Vi branching fractions for the four 
L = 1 states Q. The impact of the branching frac- 
tion uncertainties is assessed by recalculating the B —> 
D(*> (t~/£~)V PDFs and adopting the variation of the fit 
results from the ensemble of PDFs as the uncertainty. 
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FIG. 8. (Color online). Comparison of the m^iss an d \p} | distributions of the D^*'£ samples (data points) with the projections of 
the results of the isospin-unconstrained fit (stacked colored distributions). The region above the dashed line of the background 
component corresponds to BB background and the region below corresponds to continuum. The peak at m^ iss = in 
the background component is due to charge cross-feed events. The |p|| distributions show the signal-enriched region with 
m miss > 1 GeV 2 , thus excluding most of the normalization events in these samples. 



B — > D**{t /£ )v branching fractions: As noted 
above, the sharp peak in the w^iss distribution of the 
D(*'7T°£ samples constrains contributions from B — > 
D^ixlv decays. Events with additional unreconstructed 
particles contribute to the tail of the w^iss distribution 
and, thus, are more difficult to separate from other back- 
grounds and signal events. This is the case for B — » 
D**t~v t decays, which are combined with B — > D**£~vg 
decays in the D**(£/t)v PDFs with the relative propor- 
tion 1Z(D**) PS = 0.18. This value has been derived 
from the ratio of the available phase space. The same 
estimate applied to B -> D^*'£~V£ decays results in 
1l(D)p S = 0.279 and 1Z(D*) PS = 0.251, values that are 
58% and 32% smaller than the measured values. Tak- 
ing this comparison as guidance for the error on 1Z(D**), 
we increase 1Z(D**) by 50%, recalculate the D**(£/t)p 
PDFs, and repeat the fit. As a result, the values of 11(D) 
and 1Z(D*) decrease by 1.8% and 1.7%, respectively. The 
impact is relatively small, because B —> D**t~V t con- 



tributions are small with respect to signal decays, which 
have much higher reconstruction efficiencies. 

Unmeasured B — s- £>**(—!• D^n-K)ivi decays: To as- 
sess the impact of other potential B — > D**l~T>i contri- 
butions, we modify the standard fit by adding an addi- 
tional component. Out of the four contributions listed 
in Table ED the three-body decays of the D** states 
with L = 1 give the best agreement in the fits to the 
£)(*) 7r 0f samples. For this decay chain, the distri- 
bution has a long tail due to an additional undetected 
pion. This could account for some of the observed excess 
at 1 < m^ niss < 2 GeV 2 in Fig. GO We assign the observed 
change in 1Z(D^*') as a systematic uncertainty. 



2. Cross-feed Constraints 

MC statistics: Constraints on the efficiency ratios 
that link contributions from the same source are taken 
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FIG. 9. (Color online). Comparison of the m^iss and |pj| distributions of the D^n°£ samples (data points) with the projections 
of the results of the isospin-unconstrained fit (stacked colored distributions) . The region above the dashed line of the background 
component corresponds to BB background and the region below corresponds to continuum. 



from MC simulation. The impact of their statistical un- 
certainty is assessed by varying the simulated event yields 
assuming Poisson errors. 

The ratios /d**- We assess the uncertainty on /d**, 
the constraints linking the D**{t/r)v yields in the flWf 
and D^tt°£ samples, by estimating the relative efficien- 
cies of the selection criteria that differ in the two samples. 
The main differences in the selection of these samples are 
due to differences in the D^*H and D^tt°£ BDTs. 

In the D^*'l samples, we observed that differences be- 
tween data and simulation cause a 5%-10% underesti- 
mation of the continuum and BB backgrounds after the 
BDT requirements are applied. Since the D**(£/t)v con- 
tributions have similar S extra distributions, and these dis- 
tributions are the key inputs to the BDTs, we applied 
the same 5%-10% corrections to these contributions. We 
conservatively assign 100% of this correction as the sys- 
tematic uncertainty on the D**(£/t)v efficiency in the 
D^*H samples. 

Since B — > D** (r~ / 'i~)v decays are difficult to isolate 
in samples other than the D^t: £ control samples, we 



estimate the uncertainty on the D**(£/t)v efficiency due 
to the D^* s> tt°£ BDT selection by relying on the observed 
data-MC difference of the BDT selection efficiency for 
the D^tv sample. We assign the full 8.5% overestimate 
of the D(*)£v contribution as the systematic uncertainty 
on the D**(£/t)v efficiency in the D^tt°£ samples. 

The /d»« constraints also depend on the relative 
branching fractions of the four B — > D**£~~Dt decays that 
arc combined in the D**(£/t)v contributions. We esti- 
mate their impact on /d** from the branching fraction 
variations observed in the evaluation of the PDF uncer- 
tainty. The largest standard deviation for the four /jj** 
distributions is 1.8%. 

By adding the uncertainties on /d«« described above in 
quadrature, we obtain total uncertainties of 13.2% for the 
D samples, and 10.0% for the D* samples. Given that 
there are similarities between the BDT selections applied 
to the D and D* samples, we adopt a 50% correlation be- 
tween their uncertainties. With these uncertainties and 
correlations, we derive the total impact on the results, 
5.0% for TZ(D) and 2.0% for TZ(D*). 
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TABLE V. Systematic uncertainties and correlations on TZ(D^*') for the isospin-unconstrained (columns 1-4 and 7-8) and 
isospin-constrained (columns 5-6 and 9) fits. The total uncertainties and correlations are calculated based on Eg. 1271 







Fractional uncertainty (%) 




Correlation 




Source of uncertainty 


K(D°) 




TZ{D+) TZ{D*+) 


11(D) K{D*) 


D°/D*° D~ 


773*+ 


D/D* 


Additive uncertainties 




















PDFs 




















MC statistics 


6.5 


2.9 


5.7 


2.7 


4.4 


2.0 


-0.70 


-0.34 


-0.56 


B -»■ D { *\r~/r)v FFs 


0.3 


0.2 


0.2 


0.1 


0.2 


0.2 


-0.52 


-0.13 


-0.35 


D -> DW(7r /V ) 


0.7 


0.5 


0.7 


0.5 


0.7 


0.5 


0.22 


0.40 


0.53 


B(B -> D**£~V e ) 


1.0 


0.1 


1.0 


0.4 


0.8 


0.3 


-0.63 


-0.68 


-0.58 


B(B D**t~V t ) 


1.2 


2.0 


2.1 


1.6 


1.8 


1.7 


1.00 


1.00 


1.00 


D — > D { 'tt-k 


2.1 


2.6 


2.1 


2.6 


2.1 


2.6 


0.22 


0.40 


0.53 


Cross-feed constraints 




















MC statistics 


2.6 


0.9 


2.1 


0.9 


2.4 


1.5 


0.02 


-0.02 


—0.16 


/£>** 


o.z 


2.6 


5.3 


1.8 


5.0 


2.0 


0.22 


0.40 


U.Oo 


Feed- up / feed- down 


1 Q 


0.5 


1.6 


0.2 


1.3 


0.4 


0.29 


0.51 


n a 7 


Isospin constraints 




— 






1.2 


0.3 


— 


— 


— U.DU 


Fixed backgrounds 




















iviu statistics 


1 Q 

4.o 


2.3 


4.3 


1.8 


3.1 


1.5 


-0.48 


-0.05 


— U.oU 


xjiuciciicy collections 


1 8 


3.0 


4.5 


2.3 


3.9 


2.3 


-0.53 


0.20 


u. .so 


Multiplicative uncertainties 




















MC statistics 


2.3 


1.4 


3.0 


2.2 


1.8 


1.2 


0.00 


0.00 


0.00 


B -> D^ij' ir)v FFs 


1.6 


4 


1.6 


0.3 


1.6 


0.4 


00 


00 
u.uu 


0.00 


Lepton PID 


0.0 


0.6 


0.6 


0.5 


0.6 


0.6 


1.00 


1.00 


1.00 


tt'Vtt* from ZT £>tt 


0.1 


0.1 


0.0 


0.0 


0.1 


0.1 


1.00 


1.00 


1.00 


Detection/Reconstruction 


n 7 
U. / 


0.7 


0.7 


0.7 


0.7 


0.7 


1.00 


1.00 


1 on 
1.00 


£?(t~ -> l~viV T ) 


0.2 


0.2 


0.2 


0.2 


0.2 


0.2 


1.00 


1.00 


1.00 


Total syst. uncertainty 


12.2 


6.7 


11.4 


6.0 


9.6 


5.5 


-0.21 


0.10 


0.05 


Total stat. uncertainty 


19.2 


9.8 


18.0 


11.0 


13.1 


7.1 


-0.59 


-0.23 


-0.45 


Total uncertainty 


22.7 


11.9 


21.3 


12.5 


16.2 


9.0 


-0.48 


-0.15 


-0.27 



TABLE VI. Additional B -> D**£~V e decays and the MC 
model implemented for their decays. The fourth decay mode 
refers to three-body decay of the four L — 1 D" states. 



Decay 




Decay model 


Non-resonant B - 




Goity-Roberts [38] 


Non-resonant B - 




Phase Space 


B -> D M rj£jy e 




Phase Space 


B -> £>"(-> D ( * 




ISGW2 [31] 



Feed-down constraints: The feed-down constraints of 
the signal yields are corrected as part of the iteration of 
the fit. The uncertainties on these corrections are given 
by the statistical uncertainty on the ratios of the fitted 
D*£v D*l and D*£v Dl yields. They arc 2.4% and 
4.4% on the D*°tv and D* + tv feed-down constraints, 
respectively. 



Feed-up constraints: We estimate the uncertainty on 
the Dtp and D£v feed-up constraints as 100% of the 
corrections on the feed-down constraints. This results in 
6.8% on the D Q (£/r)v feed-up and 9.9% on the D+(£/t)v 
feed-up. These two effects combined lead to an uncer- 
tainty of 1.3% on %{D) and 0.4% on K(D*). 



Isospin constraints: In the isospin-constrained fit, we 
employ five additional constraints to link the signal and 
normalization yields of the samples corresponding to B~ 
and B decays. Since we reweight these contributions 
with the q 2 < 4 GeV 2 control sample, the uncertainty 
on the isospin constraints is given by the statistical un- 
certainty on the ratios of the q 2 < 4 GeV 2 yields. This 
uncertainty is 3.4% in the D£ samples and 3.6% in the 
D*£ samples. This translates into uncertainties of 1.2% 
on K(D) and 0.3% on 1Z(D*). 



22 



TABLE VII. Continuum and other BB background yields; 
the first uncertainty is due to MC statistics, the second to 
efficiency corrections, and a refers to the total uncertainty. 



Sample 


Continuum 


*{%) 


BB a 


(%) 


D°£ 


355 ± 13 ± 12 


4.9 


330 ± 6 ± 17 


5.3 


D'°£ 


132 ± 8 ± 6 


7.6 


188 ± 4 ± 10 


5.9 


D+£ 


157 ± 9 ± 6 


6.9 


191 ± 5 ± 9 


5.5 


D*+£ 


12 ± 3 ± 1 


23.6 


72 ± 3 ± 4 


6.9 



3. Fixed Background Contributions 

MC statistics: The yields of the continuum, BB, and 
cross-feed backgrounds are fixed in the fit. The uncer- 
tainty due to the limited size of the MC samples is esti- 
mated generating Poisson variations of these yields, and 
repeating the fit with each set of values. A significant 
part of this uncertainty is due to the continuum yields, 
since the size of simulated continuum sample is equiva- 
lent to only twice the data sample, 

Efficiency corrections: To account for the correla- 
tions among the various corrections applied to the con- 
tinuum and BB backgrounds, we follow this multi-step 
procedure: 

• We vary the continuum corrections within their sta- 
tistical uncertainties of 3%-9% , given by the num- 
ber of events in the off-peak data control samples. 

• The branching fractions of the most abundant de- 
cays in the BB background are varied within their 
uncertainties [l2| . 

• The BB correction is reestimatcd in the high .©extra 
control sample, and varied within the statistical un- 
certainty of 1.9%. 

• The BDT bias corrections are reestimated in the 
toes sideband, and varied within their statistical 
uncertainties, 2.1% in the D£ samples and 3.6% in 
the D*£ samples. 

• The BB background PDFs are recalculated. 

• The fit is repeated for each set of PDF and yield 
variations. 

Table EH] shows the size of the continuum and BB 
backgrounds and their uncertainties due to the limited 
size of the MC samples and the various corrections im- 
plemented by comparisons with control samples. 

B. Multiplicative Uncertainties 

MC statistics: The relative efficiency £ S ig/£norm is es- 
timated as the ratio of expected yields, so the limited size 
of the MC samples contributes to its uncertainty. We es- 
timate it assuming Poisson errors on the MC yields. 




PlI + (GeV) P7r0 (GeV) 

FIG. 10. (Color online). Pion momentum in the laboratory 
from B -»■ D* + £v and B -¥ D* + tv decays: (a) D* + -> D°tt + , 
and (b) D* + — ► D + n° decays. Histograms are normalized to 
1000 entries. 

Form factors for B -> D ( *'(r/r)i/: The q 2 > 
4 GeV 2 requirement introduces some dependence on the 
FF parameterization. This uncertainty is assessed based 
on the effect of the FF variations calculated for the un- 
certainty on the PDFs. 

7r°/7r ± from D* — > Dir: There is a significant 
momentum-dependent uncertainty on the reconstruction 
efficiency of soft pions originating from D* — > Dir decays. 
However, the momentum spectra of soft pions in signal 
and normalization decays are rather similar, see Fig. 1101 
As a result, the uncertainty on 1Z(D^) is less than 0.1%. 

Detection and Reconstruction: Given that signal and 
normalization decays are reconstructed by the same par- 
ticles in the final state, many of the uncertainties that 
impact their efficiencies cancel in the ratios e s ig/£norm- 
Uncertainties due to final-state radiation, soft-pion re- 
construction, and others related to the detector perfor- 
mance contribute less than 1%. Similarly, the tagging ef- 
ficiency for events with signal and normalization decays 
show only very small differences. 

r~ — > l~vgy r branching fraction: We use the world 
averages B(t~ -> e~V e v T ) = (17.83 ± 0.04)% and 
B(t~ -> ii-VaVr) = (17.41 ± 0.04)% [13 . 

C. Correlations 

Even though several of the uncertainties listed in Ta- 
ble |V] have the same source, their impact on TZ(D^) is 
largely uncorrelated, i.e., the correlation between uncer- 
tainties in different rows of Tabic [V] is negligible. How- 
ever, the correlation between the uncertainties on TZ(D) 
and 1Z(D*) (different columns) is significant, and impor- 
tant for the comparison of these measurements with the- 
oretical predictions. 

For most of the additive systematic uncertainties, we 
estimate the correlation from the two-dimensional 1Z(D)- 
1Z(D*) distribution resulting from the fit variations. This 
is not possible for the D** -> D^n°/7r ± and D** -> 
D^rnr uncertainties. These uncertainties affect the size 
of the D**(£/t)i> background in the D^*H samples in the 
same way that as fw does. Thus, we derive their corre- 
lations from the /r>*» correlations. Since the signal and 
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D**tv PDFs are very similar, we assign a 100% correla- 
tion on B(B -> D**t-V t ). 

The multiplicative uncertainties on the efficiency due 
to the MC statistics are uncorrelated. The FFs for 
B — > Di~vg and B — > D*l~Vi decays are measured 
separately, so their uncertainties are also not correlated. 
The uncertainty on B(t~ — > l~viv T ) affects all channels 
equally. We assume that the remaining small uncertain- 
ties on the efficiencies due to detector effects are 100% 
correlated as well. 

The uncertainties and their correlations are listed in 
Table |Vj We combine these correlations pi and the un- 
certainties by adding their covariance matrices, 



E 



*2 



PiOi(T i a. 



°tot PtotCtotCtot 
PtotCtotOtot °tot i 



(27) 

Here, Oi and a* refer to the uncertainties on 1Z(D) and 
JZ(D*), respectively. 



VIII. STABILITY CHECKS AND KINEMATIC 
DISTRIBUTIONS 

A. Stability tests 

We have checked the stability of the fit results for dif- 
ferent data subsamples and different levels of background 
suppression. 

To look for possible dependence of the results on the 
data taking periods, we divide the data sample into four 
periods corresponding to approximately equal luminosity, 
and fit each sample separately. The results are presented 
in Fig. [TTJ The eight measurements each for 7Z(D) and 
TZ(D*), separately for B + and B°, are compared to the 
isospin-constraincd fit results obtained from the complete 
data sample. Based on the values of x 2 for 7 degrees of 
freedom, we conclude that the results of these fits are sta- 
tistically consistent with the fit to the whole data sample. 

A similar test is performed for two samples identified 
by the final state lepton, an electron or a muon. This 
test includes the uncertainties on the background cor- 
rections that affect the electron and muon samples dif- 
ferently. These uncertainties are statistically dominated 
and, thus, independent for both samples. The results are 
presented in the bottom panels of Fig. [TT] The \ 2 tests 
confirm the stability of these measurements within the 
uncertainties. 

To assess the sensitivity of the fit results on the purity 
of the data sample and the BDT selection, we perform 
fits for samples selected with different BDT requirements. 
We identify each sample by the relative number of events 
in the signal region (m 2 niss > 1 GeV 2 ) with respect to the 
nominal sample, which is labeled as the 100% sample. 
The ratio of the number of fitted signal events S to the 
number of background events B varies from S/ B = 1.27 
in the 30% sample, to S/B = 0.27 in the 300% sample, 
while the backgrounds increase by a factor of 18. The 
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FIG. 11. (Color online). Measurements of K(D) and TZ(D") 
for different data subsamples. Top: for four run periods 
with statistical uncertainties only. Bottom: for electrons and 
muons with statistical and uncorrelated systematic uncertain- 
ties. The vertical bands labeled "SM" and "All data" mark 
the SM predictions and the results of the fits to the whole 
data sample, respectively. 
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FIG. 12. (Color online). Measurements of 1Z(D) and 
1Z(D*) for different BDT requirements, impacting the sig- 
nal/background ratio. The horizontal bands mark the 1Z(D) 
and 1Z(D* ) results for the isospin-constrained fit to the nom- 
inal (100%) sample. The data points represent the results of 
the fits for B + and B° mesons with their statistical uncer- 
tainties. 



BDT bias correction and the PDFs are recalculated for 
each sample. Figure H"2"1shows the results of fits to the dif- 
ferent samples with tighter and looser BDT requirements. 
We take into account the large correlations between these 
nested samples and conclude that the results are stable 
for the very large variations of the BDT requirements. 



B. Gaussian Uncertainties 

For a maximum likelihood fit with Gaussian uncer- 
tainties, the logarithm of the likelihood is described by 
the parabola P(Y) = (Y - Y Rt ) 2 /2al t , where Y &t is the 
fitted yield and am is the uncertainty on Yfi t . Figure 
[TBI compares the likelihood scan of the signal yields for 
the isospin-constrained fit with the parabola that results 
from the fitted yields, presented in Table IVIIII There is 
a slight asymmetry in the likelihood function, but good 
agreement overall. Thus, we conclude that the statisti- 
cal uncertainty on 7Z(D) and 7Z(D*) may be considered 
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Log(Likelihood) 

■ Parabola 




(N b 



1.5 GcV region. The data and simulation agree well 
within the limited statistics. 



IX. RESULTS 
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FIG. 13. (Color online). Likelihood scan for the two signal 
yields compared to a parabola. The dashed lines indicate the 
number of standard deviations (n a ) away from the fit result. 




0.45 

11(D) 



0.32 0.33 0.34 0.35 
K{D*) 



FIG. 14. (Color online). Histograms: TZ(D M ) distributions 
resulting from 1000 variations of . Solid curves: Gaussian 
fits to the TZ(D^) distributions. 



Gaussian. 

Figure [Til shows the effect on 1Z(D) and 1Z(D*) from 
variations on fw, the largest source of systematic un- 
certainty. The distributions are well described by a Gaus- 
sian function. This is also the case for the other major 
sources of systematic uncertainty. 



C. Kinematic Distributions 

We further study the results of the fit by comparing the 
kinematic distributions of data events with the SM ex- 
pectations. Specifically, we focus on the signal-enriched 
region with m^ liss > 1.5 GcV 2 and scale each component 
in the simulation by the results of the fits. To compare 
the data and MC distributions we calculate a \ 2 P er de- 
gree of freedom which only includes the statistical un- 
certainty of bins with 8 or more events. The number of 
degrees of freedom is given by the number of bins minus 
the number of fitted signal yields. 

Figure [TS] shows the i? e xtra distribution of events in the 
D(*>£ samples. This variable is key in the BDT selection 
and overall background suppression. There is a clear en- 
hancement of signal events at E ext m = in all four D^*H 
samples. The background contributions, which are sig- 
nificantly more uniform in -E xtra than those of signal, 
appear to be well reproduced. We conclude that the sim- 
ulation agrees well with the data distribution. 

Figure [TO] also shows clear signal enhancements in the 
toes an d \p*^\ distributions of events in the > 



A. Comparison with SM expectations 



Table IVIIII shows the results of the measurement of 
TZ-(D) and 1Z(D*) extracted from the fit without and with 
isospin constraints linking B + and B° decays. 

The B — > D^t~V t branching fractions are calculated 
from the measured values of TZ(D^), 

B(B -> D^t-Vt) = K(DW)xB(B -> D^rv,,). (28) 

For B~ , we use the average branching fractions measured 
by BABAR HIEl, 



B(B~ -> D°e-V e ) = (2.32 ± 0.03 ± 0.08)%, 
B(B~ -> D*°e-V e ) = (5.31 ± 0.02 ± 0.19)%, 

and for B°, the corresponding branching fractions related 
by isospin. 

We estimate the statistical significance of the measured 
signal branching fractions as £ s t&t = V /2A(EZJ, where 
A(ln£) is the increase in log-likelihood for the nominal 
fit relative to the no-signal hypothesis. The total signifi- 
cance U t t is determined as 



Cstat 



(29) 



In this expression, the statistical significance is scaled by 
the sum of the statistical uncertainty <7 8 t a t and the addi- 
tive systematic uncertainty a asys . The significance of the 
B — > Dt~T> t signal is 6.8ct, the first such measurement 
exceeding 5a. 

We compare the measured TZ(D^) to the calculations 
based on the SM, 

n(D) cxp = 0.440 ± 0.072 K(D*) cxp = 0.332 ± 0.030, 
K(D) SM = 0.297 ± 0.017 K(D*) SM = 0.252 ± 0.003, 

and observe an excess over the SM predictions for 1Z(D) 
and 1Z(D*) of 2. Oct and 2.7a, respectively. We combine 
these two measurements in the following way 



X 2 = (A,A*W w -p" 

V PC 0X p 



'th 




(30) 

where A« = ft(L>W) cxp - TZ{D^) th , and p is 
the total correlation between the two measurements, 
p(K(D),K(D*)) = -0.27. Since the total uncertainty is 
dominated by the experimental uncertainty, the expres- 
sion in Eq.[3D]is expected to be distributed as a x 2 distri- 
bution for two degrees of freedom. Figure [T7] shows this 
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FIG. 15. (Color online). E ext ra distributions for events with m^ iss > 1.5 GeV 2 scaled to the results of the isospin-unconstrained 
(first two columns) and isospin-contrained (last column) fits. The region above the dashed line of the background component 
corresponds to BB background and the region below corresponds to continuum. In the third column, the B and B + samples 
are combined, and the normalization and background events are subtracted. 
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FIG. 16. (Color online). mEs distributions before (left) and after (center) subtraction of normalization of background events, 
and lepton momentum distributions after this subtraction (right) for events with > 1.5 GeV 2 scaled to the results of the 

isospin-constrained fit. The B° and B + samples are combined. See Fig. [15] for a legend. 

TABLE VIII. Results of the isospin-unconstrained (top four rows) and isospin-constrained fits (last two rows). The columns 
show the signal and normalization yields, the ratio of their efficiencies, 1Z(D^*'), the signal branching fractions, and S s ta.t 
and I7tot, the statistical and total significances of the measured signal yields. Where two uncertainties are given, the first is 
statistical and the second is systematic. The second and third uncertainties on the branching fractions B(B — s- D'*V~77 T ) 
correspond to the systematic uncertainties due to IZ(D^) and B(B — > D^** 1 t~T>i), respectively. The stated branching fractions 
for the isospin-constrained fit refer to B~ decays. 



Decay 






-/Vnorm 


£sig/ £norm 


TZ(D M ) 


B(B D^-tv) (%) 


•^stat 




B~^f 


D°T'V r 


314 ± 60 


1995 ± 55 


0.367 ± 0.011 


0.429 ± 0.082 ± 0.052 


0.99 ± 0.19 ± 0.12 ± 0.04 


5.5 


4.7 


B~ — > 


D*°T~V T 


639 ± 62 


8766 ± 104 


0.227 ± 0.004 


0.322 ± 0.032 ± 0.022 


1.71 ± 0.17 ± 0.11 ± 0.06 


11.3 


9.4 


B° -¥ 


D + t~V t 


177 ± 31 


986 ± 35 


0.384 ± 0.014 


0.469 ± 0.084 ± 0.053 


1.01 ± 0.18 ± 0.11 ± 0.04 


6.1 


5.2 


B° —y 


D* + t~V t 


245 ± 27 


3186 ± 61 


0.217 ± 0.005 


0.355 ± 0.039 ± 0.021 


1.74 ± 0.19 ± 0.10 ± 0.06 


11.6 


10.4 


B -> 


Dt~V t 


489 ± 63 


2981 ± 65 


0.372 ± 0.010 


0.440 ± 0.058 ± 0.042 


1.02 ± 0.13 ± 0.10 ± 0.04 


8.4 


6.8 
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D*tV t 


888 ± 63 


11953 ± 122 


0.224 ± 0.004 


0.332 ± 0.024 ± 0.018 


1.76 ± 0.13 ± 0.10 ± 0.06 
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FIG. 17. (Color online). Representation of x 2 (Eq. [30]) in 
the TZ(D)-1Z(D*) plane. The white cross corresponds to the 
measured 72.(13'*'), and the black cross to the SM predictions. 
The shaded bands represent one standard deviation each. 



distribution in the 1Z(D)-1Z(D*) plane. The contours are 
ellipses slightly rotated with respect to the TZ(D)-1Z(D*) 
axes, due to the non-zero correlation. 

For the assumption that n(D^) th = H{D^) S u, we 
obtain \ 2 = 14.6, which corresponds to a probability 
of 6.9 x 10 -4 . This means that the possibility that the 
measured TZ{D) and TZ(D*) both agree with the SM pre- 
dictions is excluded at the 3.4er level [42[. Recent calcu- 
lations @, HS, S3 have resulted in values of 1Z(D) S m 
that slightly exceed our estimate. For the largest of those 
values, the significance of the observed excess decreases 
to 3.2a. 



B. Search for a charged Higgs 

To examine whether the excess in can be ex- 

plained by contributions from a charged Higgs boson in 
the type II 2HDM, we study the dependence of the fit 
results on tan/3/m#+. 

For 20 values of tan/3/m#+, equally spaced in the 
[0.05, 1.00] GeV -1 range, we recalculate the eight signal 
PDFs, accounting for the charged Higgs contributions as 
described in Sec. [H] Figure [TS1 shows the m^ iss and |p*| 
projections of the D°tv => D°£ PDF for four values of 
tan/3/m#+. The impact of charged Higgs contributions 
on the m^igg distribution mirrors those in the q 2 distri- 




gM .... tan/3/m ff + =0.5GcV 

tan/3/m H+ = 0.3GeV _1 .... tan/3/m H+ =!GeV _1 



FIG. 18. (Color online). m^ iss and projections of the 
D°rv =>• D a £ PDF for various values of tan/3 /m H + . 
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FIG. 19. (Color online). Left: Variation of the B — > Dr~V T 
(top) and B —¥ D*t~V t (bottom) efficiency in the 2HDM 
with respect to the SM efficiency. The band indicates the 
increase on statistical uncertainty with respect to the SM 
value. Right: Variation of the fitted B — > Dr~V T (top) and 
B — > D*t~V t (bottom) yields as a function of tan/3/m w +. 
The band indicates the statistical uncertainty of the fit. 

bution, see Fig. [3J because of the relation 

m miss = (Pe+e- ~ PB tB . g ~ P D l>) ~ Pi) = (? ~ Ptf , 

The changes in the |p* | distribution are due to the change 
in the r polarization. 

We recalculate the value of the efficiency ratio 
e s ig/e norm as a function of t&n(3/m H + (see Fig. I19p . 
The efficiency increases up to 8% for large values of 
tan/3/m#+, and, as we noted earlier, its uncertainty in- 
creases due to the larger dispersion of the weights in the 
2HDM reweighting. 

The variation of the fitted signal yields as a function 
of tan/3/m£f+ is also shown in Fig. 1191 The sharp drop in 
the B — » Dt~~D t yield at tan/3/m ff + ss 0.4 GeV -1 is due 
to the large shift in the m^ iss distribution which occurs 
when the Higgs contribution begins to dominate the total 
rate. This shift is also reflected in the q 2 distribution and, 
as we will see in the next section, the data do not support 
it. The change of the B — > D*t~v t yield, mostly caused 
by the correlation with the B — > Dt~v t sample, is much 
smaller. 

Figure [20l compares the measured values oiH(D) and 
K(D*) in the context of the type II 2HDM to the theoret- 
ical predictions as a function of tan/3/m^+ . The increase 
in the uncertainty on the signal PDFs and the efficiency 
ratio as a function of tan/3/m#+ are taken into account. 
Other sources of systematic uncertainty are kept constant 
in relative terms. 

The measured values of 1Z(D) and 1Z(D*) match 
the predictions of this particular Higgs model for 
tan/?/m H + = 0.44±0.02 GeV" 1 and t&n/3/m H+ = 0.75± 
0.04 GeV -1 , respectively. However, the combination of 
K(D) and K(D*) excludes the type II 2HDM charged 
Higgs boson at 99.8% confidence level for any value of 
tan/3/m#+, as illustrated in Fig. [5TJ This calculation is 
only valid for values of m H + greater than 15 GeV [1, Q. 
The region for in H + < 15 GeV has already been excluded 
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FIG. 20. (Color online). Comparison of the results of this 
analysis (light band, blue) with predictions that include a 
charged Higgs boson of type II 2HDM (dark band, red). The 
widths of the two bands represent the uncertainties. The SM 
corresponds to tan/3/m fl + = 0. 
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FIG. 21. (Color online). Level of disagreement between this 
measurement of TZ(D^*') and the type II 2HDM predictions 
for all values in the ta.nf3-m H + parameter space. 



by B — >• X s j measurements [22J, and therefore, the type 
II 2HDM is excluded in the full tan/3-m#+ parameter 
space. 

The excess in both 1Z(D) and 1Z(D*) can be explained 
in more general charged Higgs models jUSSl- The 
fectivc Hamiltonian for a type III 2HDM is 



%eff 



4G F V c i 
V2 



S L (cP L b) (tP l » t ) + S R (cP R b) (tP l v t ) , (31) 



where Sl and Sr are independent complex parameters, 
and Pl.r = (1 T 7s)/2- This Hamiltonian describes the 
most general type of 2HDM for which m 2 H+ ^> q 2 . 
In this context, the ratios TZ(D^*') take the form 

K(D) = TZ{D) SM + A D Rc(S R + S L ) + B D \S R + S L \ 2 , 
K(D*) = K(D*) SM + A D ,Rc(S R - S L ) + B D , \S R - S L \ 2 - 

The sign difference arises because B — s- Dt~V t decays 
probe scalar operators, while B — > D*t~V t decays are 
sensitive to pseudo-scalar operators. 
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FIG. 22. (Color online). Favored regions for real values of the 
type III 2HDM parameters Sr and Sl given by the measured 
values of H(D^). The bottom two solutions are excluded by 
the measured q 2 spectra. 



The type II 2HDM corresponds to the subset of 
the type III 2HDM parameter space for which Sr = 
— mf,m T tan 2 /3/m^ + and Sl = 0. 

The TZ(D^) measurements in the type II 2HDM con- 
text correspond to values of S r ±Sl in the range [—7.4, 0]. 
Given that the amplitude impacted by NP contributions 
takes the form 

\H s (S R ±S L ;q 2 )\ « \1 + (S R ±S L ) x F(q 2 )\, (32) 

we can extend the type II results to the full type III 
parameter space by using the values of TZ(D^) ob- 
tained with H S (S R ± S L ) for H S (-S R T- S L ). Given the 
small t&nf3/m H + dependence of TZ{D*) (Fig. |20|) . this 
is a good approximation for B — > D*t~~D t decays. For 
B — > Dt~V t decays, this is also true when the decay am- 
plitude is dominated either by SM or NP contributions, 
that is, for small or large values of \S r + Sl\- The shift in 
the m^jgg and q 2 spectra, which results in the 40% drop 
on the value of 71(D) shown in Fig.J^Dl occurs in the inter- 
mediate region where SM and NP contributions are com- 
parable. In this region, H S (S R + Sl) ^ H S (—S R — Sl), 
and, as a result, the large drop in TZ(D) is somewhat 
shifted. However, given that the asymptotic values of 
TZ(D) are correctly extrapolated, TZ(D) is monotonous, 
and the measured value of 1Z(D*) is fairly constant, the 
overall picture is well described by the H S (S R ± Sl) ~ 
H S (—S R T- Sl) extrapolation. 

Figure [221 shows that for real values of S R and Sl , 
there are four regions in the type III parameter space 
that can explain the excess in both 1Z(D) and 1Z(D*). 
In addition, a range of complex values of the parameters 
are also compatible with this measurement. 



C. Study of the q 2 spectra 



As shown in Sec. Ill Bl the q 2 spectrum of B — > Dt~T> t 
decays could be significantly impacted by charged Higgs 
contributions. Figure [231 compares the q 2 distribution of 
background subtracted data, corrected for detector effi- 
ciency, with the expectations of three different scenarios. 
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FIG. 23. (Color online) Efficiency corrected q 2 distributions for B — > Dt~v t (top) and B — > D*r~v T (bottom) events with 
m miss > 1.5 GeV 2 scaled to the results of the isospin-constrained fit. The points and the shaded histograms correspond to 
the measured and expected distributions, respectively. Left: SM. Center: tan/3/m H + = 0.30 GeV -1 . Right: tan/3/m H+ = 
0.45 GeV" 1 . The B° and B + samples are combined and the normalization and background events are subtracted. The 
distributions are normalized to the number of detected events. The uncertainty on the data points includes the statistical 
uncertainties of data and simulation. The values of \ 2 are based on this uncertainty. 



Due to the subtraction of the large B — > D*t~v t feed- 
down in the D£ samples, the measured q 2 spectrum of 
B — > Dt~~U t decays depends on the signal hypothesis. 
This dependence is very small, however, because the q 2 
spectrum of B — >• D*t~v t decays is largely independent 
of tan/3 /m H +. 

The measured q 2 spectra agree with the SM expec- 
tations within the statistical uncertainties. For B —> 
Dt~~D t decays, there might be a small shift to lower val- 
ues, which is indicated by the increase in the p value for 
tan/?/m H+ = 0.30 GeV" 1 . As we showed in Sec. ILTBl 
the average q 2 for tan/3/m#+ = 0.30 GeV" 1 shifts to 
lower values because the charged Higgs contribution to 
B — > Dt~T> t decays, which always proceeds via an S- 
wave, interferes destructively with the SM S-wavc. As a 
result, the decay proceeds via an almost pure P-wave and 
is suppressed at large q 2 by a factor of p 2 D , thus improv- 
ing the agreement with data. The negative interference 
suppresses the expected value of 7Z(D) as well, however, 
so the region with small tan/3/m^+ is excluded by the 
measured 1Z(D). 

The two favored regions in Fig.[22]with Sr+Sl ~ —1.5 
correspond to ta,n/3/m H + — 0.45 GeV" 1 for B — s- Dt~V t 
decays. However, as we saw in Fig. [31 the charged Higgs 
contributions dominate B — > Dt~T> t decays for values 
of tan/3/?7i#+ > 0.4 GeV" 1 and the q 2 spectrum shifts 
significantly to larger values. The data do not appear to 
support this expected shift to larger values of q 2 . 

To quantify the disagreement between the measured 
and expected q 2 spectra, we conservatively estimate the 
systematic uncertainties that impact the distributions 
shown in Fig. [231 (Appendix) . Within these uncertainties, 
we find the variation that minimizes the \ 2 value of those 
distributions. Table HXl shows that, as expected, the con- 



TABLE IX. Maximum p value for the q distributions in 
Fig. [23] corresponding to the variations due to the system- 
atic uncertainties. 







B 7Jr"77 T B 


->• D*T V T 


SM 




83.1% 


98.8% 


tan/3 /m H + 


= 0.30 GeV" 1 


95.7% 


98.9% 


ta,nf3/m H + 


= 0.45 GeV" 1 


0.4% 


97.9% 



servative uncertainties give rise to large p values in most 
cases. However, the p value is only 0.4% for B —> Dt"77 t 
decays and tan/3/m ff + = 0.45 GeV" 1 . Given that this 
value of tan/3 /m H + corresponds to Sr + Sl ~ —1.5, we 
exclude the two solutions at the bottom of Fig. [22] with 
a significance of at least 2. 9c. 

The other two solutions corresponding to Sr+Sl ~ 0.4 
do not impact the q 2 distributions of B — > Dt~V t to the 
same large degree, and, thus, we cannot exclude them 
with the current level of uncertainty. However, these so- 
lutions also shift the q 2 spectra to larger values due to the 
(S-wavc contributions from the charged Higgs boson, so 
the agreement with the measured spectra is worse than 
in the case of the SM. This is also true for any other 
solutions corresponding to complex values of Sr and Sl- 

On the other hand, contributions to B — > Dt~V t de- 
cays proceeding via P-wave tend to shift the expected 
q 2 spectra to lower values. Thus, NP processes with 
spin 1 could simultaneously explain the excess in 
[441 145| and improve the agreement with the measured q 2 
distributions. 
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TABLE X. Previous measurements of 7£(.D W ) 



Measurement 












K{D*) 






Belle 2007 [13] 










0.44 


± 0.08 


± 


0.08 


BABAR 2008 [14] 


0.42 


± 0.12 


± 


0.05 


0.30 


± 0.06 


± 


0.02 


Belle 2009 [15] 


0.59 


± 0.14 


± 


0.08 


0.47 


± 0.08 


± 


0.06 


Belle 2010 [16] 


0.34 


± 0.10 


± 


0.06 


0.43 


± 0.06 


± 


0.06 



Belle 2007 
BABAR 2008 
Belle 2009 
Belle 2010 
BABAR 2012 
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FIG. 24. (Color online). Comparison of the previous mea- 
surements of 1Z{D^) with statistical and total uncertainties 
(Table |X|| with this measurement (BABAR 2012). The verti- 
cal bands represent the average of the previous measurements 
(light shading) and SM predictions (dark shading), separately 
for 1Z(D) and 1Z(D*). The widths of the bands represents the 
uncertainties. 



by 2.0cr and 2.7er, respectively. The combined signifi- 
cance of this disagreement, including the negative corre- 
lation between 1Z(D) and 1Z(D*), is 3.4c Together with 
the measurements by the Belle Collaboration, which also 
exceed the SM expectations, this could be an indication 
of NP processes affecting B — > D^*'t~V t decays. 

These results are not compatible with a charged Higgs 
boson in the type II 2HDM, and, together with B — > X s j 
measurements, exclude this model in the full tan/3-m#+ 
parameter space. More general charged Higgs models, or 
NP contributions with nonzero spin, are compatible with 
the measurements presented here. 

An analysis of the efficiency corrected q 2 spectra of 
B —> Dt~v t and B — > D*t~T> t decays shows good agree- 
ment with the SM expectations, within the estimated un- 
certainties. The combination of the measured values of 
7£(£)(*)) anc } the q 2 spectra exclude a significant portion 
of the type III 2HDM parameter space. Charged Higgs 
contributions with small scalar terms, \Sr + Sl\ < 1.4, 
are compatible with the measured TZ(D^) and q 2 distri- 
butions, but NP contributions with spin 1 are favored by 
data. 



X. CONCLUSIONS 
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In summary, we have measured the ratios TZ(D^) = 
B(B -> D^t-V t )/B(B -> D^tTVt) based on the full 
BABAR data sample, resulting in 

11(D) = 0.440 ± 0.058 ± 0.042, 
K(D*) = 0.332 ± 0.024 ± 0.018, 

where the first uncertainty is statistical and the second is 
systematic. These results supersede the previous BABAR 
measurements (l4| . Improvements of the event selec- 
tion have increased the reconstruction efficiency of signal 
events by more than a factor of 3, and the overall statis- 
tical uncertainty has been reduced by more than a factor 
of 2. 

Table IXl shows the results of previous B —> D^t~V t 
analyses. In 2007 and 2010, the Belle collaboration mea- 
sured the absolute B — > D^t~V t branching fractions 
which we translate to TZ(D^) with B{B~ -> D°e~V i ) = 
(2.26 ± 0.11)% Q and B(B° -> D*+£-v e ) = (4.59 ± 
0.26)% For the translation of TZ{D*), we choose 

Belle's measurement of the branching fraction, instead 
of the world average, because of the current large spread 
of measured values. For Belle 2009, we average the re- 
sults for B° and B~ decays. 

The values measured in this analysis are compatible 
with those measured by the Belle Collaboration, as illus- 
trated in Fig. [Ml 

The results presented here exceed the SM predictions 
offt(£>) S M = 0.297±0.017andft(£>*) SM = 0.252±0.003 
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FIG. 25. (Color online). Assessment of the uncertainties on the g 2 distributions of background events with > 1.5 GeV 2 . 

Left: results of the isospin-constrained fit for the SM. Center: sample with 0.5 < -Ecxtra < 1.2 GeV and 5.27 < mes < 5.29 GeV. 
Right: sample satisfying the BDT requirements in the 5.20 < ro.ES < 5.26 GeV region. The data/MC plots show a fourth order 
polynomial fit and the total systematic uncertainty considered. The simulation in the control samples is normalized to the 
number of events in data. See Fig. [15] for a legend. 
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FIG. 26. (Color online). Left: g 2 distributions for the different B — > D** /£~)V contributions, all normalized to 100 events. 
Center: g 2 distributions for events with m^^s < 1.5 GeV 2 scaled to the results of the isospin-constrained fit for the SM. See 
Fig. [l5]for a legend. Right: g 2 dependence of the efficiency. The scale for the efficiency of the normalization decays is chosen so 
that the maximum value is 1. The efficiency data for the signal are adjusted so that they overlap with the data for normalization 
decays in the central part of the g 2 range. The signal efficiencies with and without the m 2 ^^ selection have the same scale. 



APPENDIX: SYSTEMATIC UNCERTAINTIES inc their sensitivity to the estimated contributions from 
ON THE q 2 SPECTRA background and normalization events. The q 2 distribu- 

tions of signal and the various backgrounds are presented 
in Fig. [25] (left). There is good agreement between the 
To assess the systematic uncertainty on the measured data and the background contributions as derived from 
q 2 distributions of B — ¥ D^*>t~V t decays, we exam- 
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the isospin-constraincd fit. To further examine the shape 
of the fixed contributions from BB and continuum back- 
ground, we show two comparisons with data control sam- 
ples: one for medium values of -E C xtra in the tties peak 
regions without the BDT requirements imposed, and the 
other for the tties sidebands with the BDT requirements. 
While the first sample shows excellent agreement over the 
full q 2 range, the smaller second sample shows some devi- 
ations at low and high q 2 . We approximate the deviation 
of the data from the simulation by a fourth order poly- 
nomial, and we adopt this difference plus the statistical 
uncertainty of each bin as the overall uncertainty of the 
BB and continuum backgrounds. We conservatively con- 
sider it uniformly distributed between the limits of the 
band shown in Fig. l25l and uncorrclated between different 
bins. 

The systematic uncertainty on the shape of the q 2 
distribution of B — > D**(t~ /i~)V decays is estimated 
by varying the relative abundance of the contributions 
shown in Fig. [251 We allow a variation of Tl(D**), the 
ratio of B ->• D**t~V t decays to B ->• D**l~Vg_ decays, 



between —20% and +50%. We also allow a contribution 
of up to 30% of B -> D**rv t decays with the D** de- 
caying into _D(*V + 7T~ . In addition, we assume a ±15% 
variation of the total B -> D**(t~ /£~)V yield. 

The q 2 spectrum of normalization decays, both well re- 
constructed and cross-feed B — > D^*>i~vn decays, is well 
described by the simulation, see Fig. [5U Given that the 
normalization decays are well understood theoretically, 
we adopt the statistical uncertainty of the simulated dis- 
tributions as the overall uncertainty of this contribution. 
Except for q 2 < 5 GeV 2 , where the rate of signal decays is 
highly suppressed, the efficiency and detector effects are 
very similar for signal and normalization. Thus, we also 
derive the overall uncertainty from the statistical uncer- 
tainty of the simulated signal q 2 distributions. 
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